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Abstract 



We consider the following stochastic partial differential equation: 

du 1 

- = -Au + KU F, 



E 



F(t,x)F(t,y) 



for x G R d in dimension d > 3, where F(t, x) is a mean zero Gaussian noise with the singular 
covariance 

5{t - s) 
\x-y\ 2 ' 

Solutions ut(dx) exist as singular measures, under suitable assumptions on the initial con- 
ditions and for sufficiently small k. We investigate various properties of the solutions using 
such tools as scaling, self-duality and moment formulae. 
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1 Introduction 



For readers who want to skip the motivation and definitions, the main results are summarized 
in Subsection |1.3| . 



1.1 Background and Motivation 



The parabolic Anderson problem is modeled by the following stochastic partial differential equa- 
tion (SPDE): 

^ = -Au + kuF. (1.1) 
dt 2 v ' 

Here u(t, x) > for t > and x G R d . The equation has various modeling interpretations (see 
Carmona and Molchanov Car94| ). The key behavior of solutions, called intermittency, is that 
they become concentrated in small regions, often called peaks, separated by large almost dead 
regions. Except when the covariance of the noise is singular at 0, the linear form of the noise 
term allows the use of the Feynman-Kac formula to study the solutions. Using this, mostly in 
the setting of discrete space with a discrete Laplacian and with a time-independent noise, there 



have been many successful descriptions of the solutions (see [GMK00] and the references there 
to work of Gartner, Molchanov, den Hollander, Konig and others.) There is less work on the 



equation with space-time noises but the memoir [ Car94 | considers the case of Gaussian noises 
with various space and time covariances. 



In addition the ergodic theory of such linear models has been independently studied. Dis- 
crete versions of the SPDE fit into the framework interacting paricle systems, under the name 
of linear systems. The reader can consult Liggett Lig85| , Chapter IX, Section 4 where, using 
the tools of duality and moments, the ergodic behavior of solutions is investigated. This work 
has been continued for lattice indexed systems of stochastic ODEs (see Cox, Fleischmann and 
Greven [ CFG96| ] and also Cox, Klenke and Perkins | CKP01| ). The basic picture is that in di- 
mensions d=l,2 and d > 3 if k is large, the dead regions get larger and larger and the solutions 
become locally extinct. Conversely in d > 3, if n is small, the diffusion is sufficient to stop the 
peaks growing and there are non-trivial steady states. 



In this paper we study a special case where the noise is white in time and has a space 
correlation that scales, namely 



E 



F(t,x)F(s,y) 



8(t - s) 
\x — y\P 



(1.2) 



The presence of slowly decaying covariances is interesting; one interpretation of the equation 



given in [Car94] is in the setting of temperature changes in fluid flow and the noise arises as 
a model for the velocities in the fluid, where it is well known that there are slowly decaying 
covariances (in both space and time). Also the equations might arise as a limit of rescaled 
models where the covariance scaling law emerges naturally. Mathematically these covariances 
are convenient since they imply a scaling relation for the solutions that allow us to convert large 
time behavior into small scale behavior at a fixed time. 
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For < p < 2 (in dimensions d > 2) there are function valued solutions with these scaling 
covariances. The Kolmogorov criterion can be used to estimate the Holder continuity of solutions 



and in Bentley | Ben99 | the Holder continuity is shown to break down as p f 2. In this paper we 
study just the case p = 2 and establish, in dimensions d > 3 and when k is small, the existence, 
and uniqueness in law, of measure valued solutions. One can imagine that the regularity of 
solutions breaks down as p j 2 but that there exists a singular, measure valued solution at p = 2 
(we do not believe the equation makes sense for the case p > 2). Note that measure valued 
solutions to an SPDE have been successfully studied in the case of Dawson- Watanabe branching 
diffusions, which can be considered as solutions to the heat equation with the noise term ^/udW, 



for a space-time white noise W (see Dawson [Daw93]). 



The special covariance \x — y\~ 2 has two singular features: the blow-up near x = y which 
causes the local clustering, so that the solutions become singular measures; and the fat tails 
at infinity which affects large time behavior (for instance we shall prove local extinction in all 
dimensions). The scaling is convenient in that it allows intuition about large time behavior to 
be transfered to results on local singularity, and vice-versa. In particular the singularity of the 
measures can be thought of as a description of the intermittency at large times. 



1.2 Definitions 



Our first task is to give a rigorous meaning to measure valued solutions of (1.1). We shall define 
solutions in terms of a martingale problem. We do not investigate the possibility of a strong 
solutions for the equation. We do however construct solutions as a chaos expansion with respect 
to any given noise. These are adapted to the same filtration as the noise and for some purposes 
provide a replacement for strong solutions. One advantage of working with martingale problems 
is that passing to the limit in approximations can be easier with this formulation. 

We now fix a suitable state space for our solutions. Throughout the paper we consider only 
dimensions d > 3. The parameter k will also be fixed to lie in the range 

0< K <^. (1.3) 

The restrictions on d and k are due to our requirement that solutions have finite second moments. 
We do not explore the possibility of solutions without second moments. 

Let «M denote the non-negative Radon measures on R d , C c the space of continuous functions 
on R rf with compact support and C% the space of functions in C c with k continuous derivatives. 
We write //(/) for the integral / f(x)fi(dx) for /jgJM and integrable /, where, unless otherwise 
indicated, the integral is over the full space H d . We consider Ai with the vague topology, that 
is the topology generated by the maps \x — > fi(f) for / £ C c . 

The class of allowable initial conditions is described in terms of the singularity of the 
measures. Define 

U\\l = I I (l + \x-y\- a )fi(dx)^dy) 
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and let H a = {fi £ M. : \\fi(dx) exp(—a\x\)\\ a < oo}. Note the spaces TC a are decreasing in a 
and increasing in a. Then define 

7~L a = [J Wq, H a+ = \j\j ftp, H a - = IJ fl ^3- 

The sets H a are Borel subset of A4. The formula for the second moments of solutions also leads, 
for each d and ft, to a distinguished choice of a. Throughout the paper we make the choice 

1/2 



d-2 



a 



d-2\ 2 



The restriction (|1.3| ) ensures a £ (0, (d — 2)/2). We shall require the initial conditions to lie in 
Ti. a +, again to guarantee the existence of second moments. 

Suppose (f2, J-~, {J~t}, P) is a filtered probability space. We call an adapted continuous A4 



valued process {ut(dx) : t > 0} a (martingale problem) solution to (IT) if it satisfies 



i. ) P(u G H a+ ) = 1, 

ii. ) {ut(dx)} satisfies the first and second moment bounds ( |1.6| ), ( |1.7| ) given below, and 

iii. ) {ut(dx)} satisfies the following martingale problem: for all / £ C% 

rt i 



zt(f)=u t (f)-u (f)- / -u s (Af)ds 

Jo 2 

is a continuous local ■Ft-martingale with quadratic variation 



<*(/)>* = ^ 



n r u s {dy)u s {dx)ds 

\x — y\ A 



(1.4) 



(1.5) 



If in addition P{uq = fj,) = 1 we say that the solution {ut(dx)} has initial condition \x. 

Let Gt(x) = {2itt)~ d l 2 exp(— \x\ 2 /2t). The moment conditions we require are that for all 
f:K d ^ [O.oo), 



E 



Mf) 



u 



G t (x - x')f(x')u (dx), 



(1.6) 



and there exists C, depending only on the dimension d and k, so that 

2 



E 



f{x)u t {dx) 



it 



C 



G t (x-x')G t (y-y')f(x')f(y') (l + 



t c 



(1.7) 

( u (dx)u (dy)dx'dy'. 



/R4d \ |x — y|"|af — y 

The construction of solutions in Section H shows that the second moment bound is quite natural. 
We believe that the moment bounds ( [1.6; ) and ( |1.7|) are implied by the martingale problem ( |1.4| ) 
and ( |1.5| ), although we do not show this. Since establishing second moment bounds is a normal 
first step to finding a solution to the martingale problem, we include these bounds as part of 
the definition of a solution. 



We finish this subsection with some simple consequences of the second moment bound. 
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Lemma 1 Suppose {ut(dx)} is a solution to with initial condition /i. Choose a so that 

fi(dx) exp(— a\x\) £ TL a . 



i) For any f £ C c and t > 0, we have 

rt r <■ f(x)f(y) 



E 



/o 



\x -y\< 



u s (dx)u s (dy)ds 



< oo 



and hence the process zt(f) defined in (E^) is a true martingale 
ii) For any < p < d — a and t > 

(1 + \x - y\~ p ) e~ a ^- a ^u t {dx)u t (dy) 



E 



< oo 



and hence ut £ TC^ a )- almost surely. 



Proof. For part i) it is sufficient to check that E[{z(f))t] < oo to ensure that zt(f) is a true 
martingale. Using the second moments (1.7) we have 



E 



f(x)f(y) 



u s (dx)u s (dy)ds 



Uo J J \x- y\- 

< eft G s (x-x')G s (y-y') J ^ JJK " ; (l+ , 
Jo jRid \x'-y'\ 2 V \x-y\ a \x'-y' 



(1.8) 



j-; n(dx)fi(dy)dx dy ds. 



We now estimate the dx' dy' integral in the above expression by using the simple bound, for 
< r < d, 



G t (x - x')G t (y - y')\x' - y'^dx 1 dy' < C(r) (|x - y\- r A t^' 2 ) . (1.9) 
For compact support / and any a > we have the bound 

J G s (x- x')f(x')dx' < C{a, f, t)e~ a ^ for all s < t, x £ K d . (1.10) 
Then, applying Holder's inequality with 1 < p < d/2 and p~ l + q~ l = 1, we have, for all s < t, 

G s (x - x')G s {y - y') \ {x y { f.} dx'dy' 
\x — y'\ z 

G s {x - x')G s {y - y')— j-^-dx'dy' 

\x — y'rP J 



< 



y'\ 2p 

G s (x - x')G s (y - y')P{x')n y ')dx'dy' 

< C(a,f,t)e- a ^- a ^ (Ix-yl^As- 1 ) 

< C{a, f, t)e- a ^- aM \x - y\~ aA2 s ' { - 2 ~ a) + /2 

< C{a,f,t)e- a ^- aM (l + \x-y\- a ) s -( 2 -«)+/ 2 . 



1/9 
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A similar calculation, using 2 + a < d, gives the bound 



G s {x - x')G s (y - y')J^M^ldx'dy' < C(a, f,t) S -^ 2 



F - y 

Now we substitute these bounds into ( |1,8| ) to obtain 

f(x)f(y) 



E 



F - v\ 



-u s (dx)u s (dy)ds 



< C(a,f,t)j s -(2-«)+/ 2 / / (l + \x-y\- a )e- a W- aM n(dx)n(dy)ds. 



which is finite since \\p(dx) exp(— a|x|)|| Q < oo. 

For part ii) use the second moment bound ([O]) to see that 



E 



(1 + \x - y\~ p ) e- alxl - aM u t (dx)u t (dy) 
< C [ G t (x - x')G t (y - y') (l + \x' - y'\- p ) e^-^ 



(1.11) 



1 + 



\x — y\ a \x' — y 



^ I p(dx)p(dy)dx dy . 



Using the bound J G t (x — x') exp(— a\x'\)dx' < C(t, a) exp(— a\x\) and ( |1.9| ), we estimate the 
dx'dy' integral in a similar manner as above. We illustrate this only on the most singular term. 
For p, q > 1 with p^ 1 + q^ 1 = 1, 



G t (x - x')G t {y - y')\x' - y'^^+^e-^'^'Ux' dy' 
G t (x - x')G t {y - y')e- ap \ x '\- a P\y'\ dx'dy') 
• (I J G t {x - x')G t {y - y')\x' - y'\-^ +a Ux'dy' 



< 



1/9 



< C(a,p,p,q)e- alxl - alvl (\x-y\- {p+a) At-( p+a ^ 2 ) (1.12) 

< C(t,a,p,p,q)e- a ^- a ^. 

provided that q(p + a) < d. Such a q > 1 can be found whenever p + a < d. Substituting this 



estimate into (1-11) gives the result. 



1.3 Main Results 



We start with a result on existence and uniqueness. 

Theorem 1 For any p £ T~t a + there exists a solution to started at p. Solutions starting 

at p G Ti. a + are unique in law. If we denote this law by then the set {Q^ ■ p S 7~La+} forms 
a Markov family of laws. 
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The existence part of Theorem |l| is proved in Section ||] and the uniqueness in Section ||. The 
next theorem, which is proved in Section ||, shows death from finite initial conditions and local 
extinction from certain infinite initial conditions. Write B(x,r) for the open ball or radius r 
centered at x. We say that a random measure uq has bounded local intensity if E[uq(B(x, 1))] 
is a bounded function of x. 

Theorem 2 Suppose {ut(dx)} is a solution to ( \1. ]\ ). 

i) Death from finite initial conditions. If P{uq G H^,) = 1 then (u t , 1) — ► almost surely 

as t — > oo. 

ii) Local extinction from infinite initial conditions. If uq has bounded local intensity and 

A C R d is a bounded set then ut(A) — > in probability as t —* oo. 

Finally, we state our main results describing the nature of the measures ut(dx). These are 
proved in Section [||. 

Theorem 3 Suppose that {ut(dx)} is a solution to \l.tp satisfying P(uq ^ 0) = 1. Fix t > 0. 
Then the following properties hold with probability one. 

i) Dimension of support. If a Borel set A supports the measure ut{dx) then the Hausdorff 

dimension of A is at least d — a. 

ii) Density of support. The closed support of ut(dx) is R d . 

iii) Singularity of solutions. The absolutely continuous part of ut(dx) is zero. 



Remarks 

1. Although Theorem || gives an almost sure result for fixed t, it leaves open the possibility 
that there are random times at which the properties fail. In Section || we shall show that 
P(u t S T~t a + for all t > 0) = 1. This implies that the weaker lower bound d — 2 — a on the 
dimension of supporting sets is valid for all times. 

2. The reader might compare the behavior described in Theorem || with that of the Dawson- 
Watanabe branching diffusion in R rf , for d > 2. This is a singular measure valued process whose 
support is two dimensional, and, if started with a finite measure of compact support, has compact 
support for all time. 

3. Many of the results go through for the boundary case K = (d — 2)/2 and for initial 
conditions in Ti a , although we have not stated results in these cases. The chaos expansion in 
Section |3| holds in both these boundary cases and the second moments are finite. Although our 



proof that the chaos expansion satisfies (|l.5[) uses k < (d — 2)/2 and /i G H a + we do not believe 



these restrictions are needed for this. However our proof of uniqueness for solutions in Section 
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H does seem to require the strict inequalities. This leaves open the possibility that there are 
solutions with a different law to that constructed via the chaos expansion. Theorems ^ and 
H will hold in the boundary cases for the solutions constructed via chaos expansion solutions. 
Parts of Theorems ^ and || also hold for all solutions, for example Propositions |2| and || use only 
the martingale problem in their proof and hold for any solution in the boundary cases. 



1.4 Tools 



We briefly introduce the main tools that we use. The first tool, simple scaling for the equation, 
is summarized in the following lemma. 

Lemma 2 Suppose that {ut(dx)} is a solution to \ ). Let a, b, c > and define 

v t (A) = au bt {cA) for Borel i4CR d 

where cA = {cx : x £ A}. Then {vt(dx)} is a solution to the equation 

dvt b . b 1 / 2 ■ . , 
- w = ^Av + K —vF biC (t,x) 

where Ff) C (t,x) is a Gaussian noise identical in law to F(t,x). 



The equation for {vt(dx)} is interpreted via a martingale problem, as in fll.ip . The easy proof 
of this lemma is omitted. 



The next tool is our equation for the second moments. The linear noise term implies that 
the solutions have closed moment equations. By this we mean that the moment densities 

H t (x 1 ,x 2 , ■ ■ ■ ,x n )dxi . . . dx n = E [ut{dxi)u t {dx 2 ) ■ ■ ■ u t (dx n )] 

satisfy an autonomous PDE. Formally assuming the solution has a smooth density ut(x), ap- 
plying Ito's formula to the product ut(x\) . . . Ut{x n ) and taking expectations suggests that H t 
satisfies 



Ot 



2 ' 



l<i<j<n I""' *JI 

Then the Feynman-Kac representation for this linear equation suggests that 



Hf\XXj • • • , Xn) F xl ^ ^ Xn 



u (X})...u (X?)exp( f 

Wo ^fri. 



l<i<j<n \ A s ^ s I 



-ds 



where E Xl ^__ jXn denotes expectation with respect to n independent d-dimensional Brownian 
motions. This formula makes sense when uq has a density, but more generally we can expect 
for solutions {u t (dx)} to ([O]) started at fi, and when /j G C c for i = 1, . . . , k, 



E 



R 2nd 



Tj4,yi,...,y„ 
0,xi , . . . 



cxp 



E 



l<j<k<n 



X. 



(i) 



X. 



(k) 



rds 



(1.13) 

n 

Yl G t (xi - yi)fi{yi)ii{dxi)dyi 



i=i 
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where Eq V ^ '"'^™ is expectation with respect to n independent d-dimensional Brownian bridges 

(X^\ . . . , X t ) started at (xj) at time zero and ending at (jji) at time t. In Section || we 
investigate the values of k for which this expectation is finite. 

The next tool is the expansion of the solution a Wiener chaos expansion, involving multiple 
integrals over the noise F(t,x). Wiener chaos expansions have been used before for linear 
equations; for example see Dawson and Salehi | Daw80| 1 or Nualart and Zakai NZ89 |. The idea 



is to start with the Green's function representation, assuming (falsely) that a function valued 
solution exists: 

u t (y) = Gtn(y) + k J J G t - 8 {y ~ z)u s (z)F(dz,ds). (1.14) 

The first term on the right hand side of this representation uses the notation Gt^{y) = J Gt(y — 
z)n(dz). The second term involves again the the non-existent density u s (z). However we can 
use the formula for ut(y) given in ( 1.14j) to substitute for the term u s (z) which appears on its 



right hand side. The reader can check that if we keep repeating this substitution, and assume 
the remainder term vanishes, we will arrive at the following formula: for a test function / G C c , 

oo 

Mf) = £ 4 n) (/./i) (i.i5) 

n=0 

where 

4 n) (f,v) =11 f(y)4 n) (y,zMdz)dy (1.16) 



and where the 1^ are defined as follows: I^°\y, z) = Gt(y — z) and for n > 1 

4: } +1 (yn + i,z) = K n f n ... f 2 [ G sl (y 1 -z)f[G Si+1 „ Si (y i+1 - yi )F(dy i ,ds i ). (1.17) 

+ Jo Jo Jo Jn nd fj[ 



In Section |3| we shall show that the stochastic integrals in (1.17) are well defined, and the series 
( 1.15D converges in I? and defines a solution. The point is that the series J2 n ^t^iVi z ) does not 



converge pointwise, but after smoothing by integrating against the initial measure and the test 
function the series does converge. The restriction fll.3| ) on k and the choice of space H. a + for 
the initial conditions is exactly what we need to ensure this L 2 convergence. For larger values 
of k it is possible that the series converges in LP for some p < 2. It is also always possible to 



consider the chaos expansion (1.15) itself as a solution, if we interpret solutions in a suitably 
weak fashion, for example as a linear functionals on Wiener space. We do not investigate either 
of these possibilities. 

The symmetry of the functions l[ n \y,z) in y and z makes it clear that a time reversal 
property should hold. This is well known for linear systems and for the parabolic Anderson 
model, and is often called self duality. Suppose that {ut(x)}, {vt(x)} are two solutions of Jl.lj) 
started from suitable absolutely continuous initial conditions uo(x)dx and vo(x)dx. We expect 
that ut(vo) has the same distribution as uo(vt). In Section Q we shall use this equality to establish 
uniqueness of solutions. 

The Feynman-Kac formula is a standard tool in analogous discrete space models. In the 
continuous space setting of the parabolic Anderson equation (01) , we shall replace the noise F 
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by a noise F that is Gaussian, white in time and with a smooth, translation invariant covariance 
T(x — y) in space. Then the Feynman-Kac representation is 



u t {x) 



E x 



G t (x-y)u (y)e- r ^ t E t % 



exp ( k I F(ds,X t _ s 



exp ( k I F(ds, X s 



(1.18) 



A proof of this representation can be found in Kunita [Kun9C] Theorem 6.2.5 and we make use 
of it in Section ||. Since our covariance blows up at the origin the exponential factor T(0) is 
infinite and the representation can only be used for approximations. 

Finally a remark on notation: we use C(t,p, . . .) for a constant whose exact value is unimpor- 
tant and may change from line to line, that may depend on the dimension d and the parameter 
k (and hence also on a), but whose dependence on other parameters will be indicated. 



2 A Brownian exponential moment 

As indicated in the introduction, the second moments of solutions {u t (dx)} to fll.lD can be 
expressed in terms of the expectation of a functional of a Brownian bridge. An upper bound 
for these expectations is a key estimate in the construction of our solutions. In this section we 
show the following bound. 



Lemma 3 For all < r\ < (d — 2) 2 /8 there exists C(rj) < oo so that for all x, y, t 



jpt,y 



where 



exp 7] 



a{rj) 



ds 

JxJ- 



d-2 



< C( V ) (1 + ■ 
d-2\ 2 

- 2rj 



t 



x \ \y\ 

1/2 



We first treat the case of Bessel processes and Bessel bridges (see Rcvuz and Yor [ RY91 1 chapter 
XI for the basic definitions). The reason for this is that the laws of two Bessel processes, of 
two suitable different dimensions, are mutually absolutely continuous and the Radon-Nikodym 
derivative involves exactly the exponential functional we wish to estimate. 

Let C[0, t] be the space of real valued continuous paths up to time t and let {Rt} be the 
canonical path variables. For d £ [2, oo) and a, b > we write E^ 1 for expectations under the 
law of the cf-dimensional Bessel process started at a and qf^ (a, b) for the transition density. We 
write E^fl t for expectations under the law of the ci-dimensional Bessel bridge starting at a and 
ending at b at time t. Suppose that Y is a non-negative random variable on the space C[0, t], 
measurable with respect to cr(R s : s <t). Lemma 4.5 of Yor [ Yor80fl , (or Revuz and Yor |RY9l| 1. 
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Chapter XI, exercise 1.22), expressed in our notation, states that the following relationship 
holds: if A, \i > then 

-Al 

= _g;(2M+2) 



,2 rt 



2 Jo R 2 a I V a 



E {2\+2) 

a - . ^ 

Now for < rj < (d — 2) 2 /8 we choose values for A, /i, Y in this identity as follows: 



(2.1) 



, d-2 



d-2 



nl/2 



2r/ 



, y = exp 



77 H 

' 2 



Note with these choices that a(rj) = A — /x, 2r] + /j 2 — A 2 = 0, and ci = 2A + 2. Applying ( |2TT| ) 
we find 



E (d) 



,2 



B (2A+2) 
-■-'a 

E {2fx+2) 
E {2^+2) 



A 2 f* ds\ [R t 



Y exp 



m 



exp 



(j? A 2 

77 H 

' 2 2 



l(i2t G c?&) 



o- a W6 a Wg| 2A,+2) (a,6)d6. 



Hence 



a,b,t 



exp 77 



(is \ 



(2^+2) / 7 \ 
0(77) ^a(T)) gt l a ' °J 



(2.2) 



/o « 2 y 

There is an exact formula for the Bessel transition density 

q { t d \a, b) = rV^-W/ 2 exp(-(a 2 + tf)/2t)I myi (ab/t) 

in terms of the (modified) Bessel functions I„ of index 1/ = (d/2) — 1. The Bessel functions l v {z) 
are continuous and strictly positive for z £ (0, 00) and satisfy the asymptotics, for ci, C2 > 0, 

I„(.z) ~ ciz^ as z J, 0, ~ C22 _1 / 2 e 2: as z f 00. 

Using these we find that 

* Y 



GXPl??/ o 



< C(n) + - for all a,b,t > 0. 



(2.3) 



We now wish to obtain a similar estimate for a Brownian bridge. Recall the skew product 
representation for a d-dimensional Brownian motion Xt, started from x 7^ 0. There is a Brownian 
motion W(t) on the sphere S d_1 , started at x/\x\ and independent of X, so that 



X t /\X t \=W[ j IAJ--V, ') . 
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We may find a constant C so that P x {W{t) £ dO) < Cd9 for all x G S d-1 and i > 1. We now 
consider the exponential moment for a <i-dimensional Brownian bridge running from x / to 
y ^ in time t. 



exp 7/ 



o 



ds 

\x7\- 



.c.O 



exp 







ds 



ds 

JxJ- 



> 1 



(2.4) 



Now we estimate the second term on the right hand side of ( |2.4| 
* ds \ _ ( ft ds 



F y,t 



exp r] 



XI 2 



< 



1 


G t (x- 


y) 




-d 


G t (x- 


y) 




-d 



\A-s 



Gt(x-y) 



1 

exp 
exp I 7] 



exp 7] 



\X S \ 2 
1 ds 

\X I 2 
K v s I 

* ds 

o W 
1 ds 

o 



> 1 



o 



* ds 

\X I 2 



G t (x-y) H'W'* 



exp 77 



ds 
\X I 2 

ds 



> l,X t G 



> l,|X t | €c%|,W 



> 1,1^1 G% 



(is 



l-fts 



Using the explicit representation for the Bessel density given above we find that 

\y\ 1 - d tf\\AM) <C{R) wh enever < R. 

Gt(x-y) ~ v f " 



Combining this with ( |2.4| ) and our estimate for the Bessel bridge we obtain the desired 

bound for (x,y,t) in any region where {|a;||y|/i < R}. 



We felt there should be a short way to treat the remaining case, but we seem to need a 
slightly complicated argument to treat the case large. Note our aim is only to find a 

constant bound for the exponential moment in this region. Define 



F K {x,y,t) 



exp 77 



* ds 
\XJ- 



A K 



< E t,y 



exp 77 



ds 



--: F(x,y,t). 



Brownian scaling implies that F(x, y, t) = F^^x, c 1//2 y, ct) for any c > 0. So we may scale 
time away and it is enough to control F(x,y, 1). We have proved above, for any R, 



F(x,y,l) <C(R, V ) 1 + 



m \y\ 



whenever \x\\y\ < R. 



(2.5) 



We first show we may reduce to the case where \x\ = \y\. Suppose that \x\\y\ > 1 and |x| > \y\. 
Define stopping times 

o x = inf{t : \X t \ < \y\}, a 2 = mf{t : \X t \\y\ < 1 - 1}, 

and let a = u\ A oi- Note that for t < o\ we have 1/|X| < l/|y| and for t < 02 we have 
1/1-^*1 < |?/|/(1 ~~ t). So we can bound the integral in F(x,y, 1) by 
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dt r dt f 1 dt 



" Jo |X t |2 J fj |X t |2 







(1-t) 2 |y| 2 / ^ |A^| 2 



Kl-|2/| 2 )+ | y |2 rl I i C 1 dt 



< 2 + 



(i-t)2 Vm 2 ) + lyl 2 ^ l^l 2 

1 



2 ' 



Conditioned on the values of a and X a , the path between t G [cr, 1] is a new Brownian bridge. 
Hence 



F(x, y, 1) < e 2r >E (F(X a , y,l-a))=E 



F[ V - 1 

1 (l-ax) 1 / 2 ' (1-^)1/2' 



By definition \X a2 \\y\/(l - a 2 ) = 1 so that F(X a2 /(l - a 2 ) 1/2 , y/(l - cri) 1/2 , 1) can be bounded 
by a constant using fl2.5| ). Also on the set {<ti < 02} we know that |X CTl ||y|/(l — o\) > 1 and 
l-^o-i | = \y\- So if we can bound the F(x,y,l) on the set diagonal case {|x| = |y|,|x||y| > 1} we 
can bound F(X ai /(l — (j\) 1 / 2 ,y / (1 — ai) 1 / 2 ,1), and in consequence also F(x,y,l). 

We now give a brief sketch to motivate the final argument. Consider the "worst case" of a 
bridge from x = Ne\ to y = —Ne\ over time one. Run both ends of the bridge until both ends 
first hit the ball of radius N/2. When TV" is large the bridge will enter the ball near x/2 and exit 
near y/2 and spend close to time 1/2 inside the ball. We may therefore approximately bound 
the exponential as 



cxp 



Using the scaling of F(x,y,t) we see that F(Ne\, —Ne\, 1) is approximately bounded by 

exp(4iV- 2 )F(iVei/2, -Nei/2, 1/2) = exp(4iV^ 2 )F(iVei/2 1 / 2 , -Ne l /2 l l 2 , 1). 

By iterating this argument we will bound F(Ne\, —Ne%, 1) for large N by values for small N 
where we know it is bounded by (2.5). 

We now give the basic iterative construction. Suppose that |x| = \y\ = R > 1 and consider 
the Brownian bridge {Xt} from x to y in time 1. Define random times 

a = inf{t : \X t \ < R/2}, r = sup{t : \X t \ < R/2} 

on the set {inf t \X t \ < R/2} = {a < r}. On the set {\X t \ > R/2, Vt G [0, 1]} we have the bound 
/o \X s \~ 2 ds < \R- 2 . On {a < r} we have the bound 

1 dt < 4R~ 2 r dt 



\X,\ 2 ~ L \XtP- 
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Conditioned on cr, r, X a , X T , the part of the path {X t : t G [<r, r]} is a new Brownian bridge. So 
we may estimate 

F(x,y,l) <exp(4r]R- 2 )(P({\X t \ > R/2, Vi G [0, 1]}) + £ (F(X CT ,X r ,r - a)l(a < r))) . (2.6) 
The same bound holds with F replaced by Fk- 

We will repeat this construction with a new Brownian bridge running from X a j{r — a) 1 / 2 
to X t /(t — a) 1 / 2 . The following lemma shows that when R is large we have usually made an 
improvement in that this bridge is closer to the origin. 

Lemma 4 There exists 7 < 1 and C3 < 00, C4 > so that, when \x\ = \y\ = R, 

P(a<T, X a -X T < 0, \x a /(T - a) l/2 \ > 7 |s|) < c 3 exp(-c 4j R), (2.7) 
and there exist C5 < 00, cq > so that, if in addition x ■ y > 0, 

P ( inf \X t \ < R/2) < c 5 exp(-c 6j R), (2.8) 



Proof. We scale the Brownian bridge by defining Xf" = Xt/R. The starting and ending 
positions x = X^, y = Xf- now satisfy |x| = \y\ = 1 and the process X^ is stopped upon hitting 
the ball of radius 1/2. However the process X^ has reduced variance. Indeed, in law we have 
the equality 

Xf- = (1 - t)x + ty+ (B t - tB^/R. 

As R — > 00 the process converges to the straight line Xt = (l — t)x + ty. For this limiting process 
the basic construction is deterministic. If x ■ y > then the straight line never gets closer to the 
origin than 2 -1 / 2 . For large R a large deviations estimate shows that deviations away from the 
straight line are exponentially unlikely and ( |2.8| ) follows. To obtain ( |2.7| ) one again considers 
the straight line Xt and maximizes X a j{r — a) 1 / 2 over those starting and ending points x,y 
for which X a ■ X T < 0. The maximum occurs, for example, when X a /(T — cr) 1 / 2 = (l/2)ei 
and X T j(r — cr) 1 ' 2 = (l/2)e2- A little trigonometry show that either X a ■ X T < or else 
X c /(t — a) 1 / 2 < 7 for some 7 G (0, 1). By taking 7 G (7, 1) a large deviations argument yields 



Applying (|2.Sj ) to the bound (|2.6|) we find, when \x\ = \y\ = R and x ■ y > 0, 

F K {x,y,l) <exp(4 V R- 2 )ll + c 5 e- C6R sup F K {x,y,l)) . (2.9) 

V W=\y\>R/2 J 



Now we wish to iterate the basic construction to define a Markov chain (x(n), y(n)) n= o,i,... 
on ^H d U {A}^ . Throughout |x(n)| = |y(n)| or x{n) = y(n) = A will hold. A is cemetery state 
from which there is no return. It will be convenient to set -F(A, A, 1) = i^(A, A, 1) = 1. We 
set x(0) = x,y(0) = y. Suppose x(n),y(n) have been defined and are not equal to A. Then we 
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repeat the basic construction described above, but started at the radius R = \x(n)\ = \y{n)\. 
We define 

x(n + 1) = XJ(t - a) 1 / 2 , y(n + 1) = X t /(t - a) 1 ' 2 on {a < r}, 
x(n + 1) = y(n + 1) = A on {\X t \ > R/2, Vt G [0, 1]}. 



We will choose a constant Ro € [1,R] shortly. Define stopping times for (x(n),y(n)) as 
follows 

Ni = inf{n : x(n) = y(n) = A}, 
N 2 = inf{n : \x{n)\ < R }, 
N3 = inf{n : x(n) ■ y(n) > 0}, 
N4 = inf{n : \x(n)\ > "f\x(n — 1)|}. 

Let N = Ni A N2 A N3 A N4. Technically we should define | A| to make these times well defined, 
but we adopt the convention that if N > k and N\ = k then N2 = N3 = N4 = 00. Note that N 
is a bounded stopping time since if N4 has not occurred then N2 < Nq, where Ra N ° < Ro. We 
now expand Fjc(x,y, 1) as in Q2.6| ) to find 

N 

F K (x,y, 1)<Y,E [l(N = n) exp ^ V (\x(0)\- 2 + ... + \x(n - 1)|" 2 )) F K {x(n),y{n), 1) . 

n=l 

On = n} we know that i?o < \x{n — 1)| < 7|x(n — 2)| < 7 2 |x(n — 3)| < ... Hence on this set 

4?/ 



exp (4j?(|a;(0)r 2 + . . . + \x(n - 1)|- 2 )) < exp (j-^ 



7 2 : 



We choose Rq large enough that this exponential is bounded by 2. This leads to the simpler 
bound 

F K (x,y,l) < 2E[F K (x(N),y(N),l)} (2.10) 

We now find various estimates for E [Fk(x(N) , y(N) , 1)] depending on the value of N. When 
N = N\ we have, by definition, 

l(JV = JVi)i^(a!(JV),y(^),l) = l. (2.11) 

When N = N2 we have |x(n)| G [Rq/2, Rq] and so we can bound 

1(N = N 2 )F K (x(N),y(N),l) < sup F(x,y,l). (2.12) 

\x\=\y\e[Ro/2,Ro] 

When N = N 3 > N 2 we have \x(N)\ = \y(N)\ > R and x(N) ■ y(N) > and we may use (|J) 
to bound FK(x(N),y(N), 1). By choosing Rq large enough this gives the bound 

1(N = N 3 > N 2 )F K (x{N),y(N),l) <2[l + -}- sup F^(x,2/,1)|. (2.13) 



16 \x\=\y\>Ro/2 



Finally when N = N4 < N2 A N3 we simply bound 



E[1(N = N 4 <N2AN 3 )F K (x(N),y(N),l)} < P(N = N 4 < N 2 V N 3 ) sup F K {x,y,l). 

\x\=\y\>R /2 
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We now claim that 



lim sup P(N = N 4 < N 2 A N 3 ) 

R^OO \ x \ = \y\>R 



0. (2.14) 

Indeed we may apply Lemma |] to see that 

P(N = N 4 = k + 1< N 2 A N 3 \(x(j),y(j)) j = 0, 1, . . . , k) < c 3 exp(-c 4 |x(A;)|)l(iV > k) 

So 



P(N = iV 4 = k < N 2 A N 3 



k=l 

oo 

< ^TE[c 3 exp(—c 4 \x(k — 1)|)1(JV > fc — 1)] 

k=l 

- N-l 

< E c 3 ^2 exp(-c 4 |x(/c)|) 

k=0 

oo 

< c 3 ^2 exp(-c 4 i?o7~ fc ) 

fe=0 

— > as i?o — > oo. 
Using the claim ( |2.14D we may choose Rq large enough that 

E[l(N = N 4 >N 3 VN 2 )F K (x(N),y(N) 1 l)]<l sup F^(x,y,l) 

B |o;|=|s/|>flo/2 



(2.15) 



Choosing i? large enoug h that all four estimates (pUP , (|1|), (p7T3|) , (gig) hold, we substitute 
them into ( 2.1C| ) to obtain 



^(x.J/.l) < 6 + 2 



< 6 + 3 



sup 

M=|2/|e[#o/2,-Ro] 
sup 

kl=|j/|6[Ro/2,i?o] 



F(x,y,l) + 



sup F K (x,y,l) 

\x\ = \y\>R /2 



1 



1 M=M>«o 



Take the supremum over x, y in {\x\ = \y\ > i?o} of the left hand side to obtain 

sup Fx(x,y, 1) < 12 + 6 sup F(x,y,l). 

M,ll/|>-Ro/2 [a:[=[v[e[iiQ/2 ) Jio] 

Letting — > oo we have bounded F(x,y, 1) on the set {|x| = |y| > i?o/2}. Together with ( |2.5| ) 
this completes the proof of the main estimate. 



Remarks 



1. The moment E 



(d) 

a,b,t 



exp(rj J"q R s 2 ds) is infinite for r) > (<i — 2) 2 /8. This follows since the 
formula ( |2.2j ) cannot be analytically extended, as a function of r/, into the region {z : Re(z) < r} 
for any r > (d — 2) 2 /8. This strongly suggests there are no solutions to (1.1) having finite second 
moments E[(ut (/)) 2 ] when k > (d — 2)/2. Similarly, the blow-up of the Brownian exponential 
moment suggests there should be no solutions to ( |1.1[) with finite second moments for any k > 
when the noise has covariance (|1.2[) with p > 2. 
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2. As indicated in Subsection 1.4 , higher moments are controlled by the Brownian expo- 
nential moments ( p.. 13 ). Using Holder's inequality we find 



pt,yi,...,y n 



exp 



l<j<k<n 



yU) y-(fc) 



rds 





( 




< n 




exp J 


l<j<k<n 







n(n — 1)k 2 



rds 



l/n(n-l) 



y(i) 



The exponential moment calculated in this section shows that this is finite when n(n — l)n 2 /2 < 
(d — 2) 2 /8. This should lead to the solutions to (1.1) having finite moments E[(ut(f)) n ] when 
k < (d — 2)(4n(n — l))" 1 / 2 . We do not think this simple Holder argument leads to the correct 
critical values for the existence of higher moments. 



3 Existence of Solutions 



In this section we give a construction of solutions to ( |l.l[) using the chaos expansion ( |l,15| ). 
However, it is hard to show from the series expansion that the resulting solution is a non- 
negative measure. For that purpose we give a second construction as a limit of less singular 
SPDEs. A comparison theorem will show that the approximating equations have solutions 
which are non-negative functions implying that the limit must also be non-negative. Finally we 



show that the two constructions yield the same process and that it is a solution of (1.1) 



We first construct a noise F with the desired covariance. Let g{x) = C7\x\~( d+2 ^ 2 . A simple 
calculation shows, for a suitable value of the constant C7, that the convolution g * g{z) = \z\~ 2 . 
Now let W be an adapted space-time white noise on n d x [0, 00) on some filtered probability 
space (O,.? 7 , {^i}, P). Define, for / : H d — ► R that is bounded, measurable and of compact 
support, 

F(t, /) = // (/ * 9)(z)W(dz, ds). (3.1) 
It is straightforward to show that F(t, f) is well defined, is a Gaussian martingale, and that 

mm 



\y 



-dy dz. 



If we write F(t, A) when / = I a then {F(t, A) : t > 0, A C H d } is a martingale measure and 
hence (see | Wal86 | Chapter 2) can be used to define a stochastic integral Jq J h(s,y)F(dy,ds) 



for suitable predictable integrands h so that 

h(s,y)F(dy,ds) = 



h(s,y)h(s,z) 



dy dz ds. 



Next we show that the expansion ( [1.15D converges. 
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Lemma 5 Suppose \i G Ti a - Then, for f G C c , the series Y^=o ^t^C/j A 4 )' defined by ( Lit ) and 
fitl.lTQ , converges in I? . Moreover 



E 



\n=0 / 



?1=0 



4 B) (/,/*) 



R 4d 



f(3/)f(a/)G t (x-2/)G t (y-T/) 



(3.2) 



j?t,x',y' 



exp 



dx'dy' fj,(dx)fj,(dy). 



Proof. We first check that the right hand side of ( |3.2| ) is finite. Using the fact that {X} — X|) / \/2 
is a Brownian bridge from x — y to x' — y' we may use Lemma to obtain 



R 4d 



fi3/)f^)G t {x-^)G t {y-t/)^ 



x.y 



exp 



o |JO-X 2 | 



fj,(dx)fj,(dy)dx'dy' 



< C(t) / /(y')/(ar')G t (x - x')G f (y - y') (1 + |x - y\- a \x' - y'\~ a ) ^{dx)^{dy)dx' dy' 
Now estimates as in Lemma |l| show this expression is finite. 



The multiple Wiener integrals of different orders are orthogonal, if they have finite second 
moments; that is if m ^ n, and if 



E 



< oo 



(3.3) 



for k = m,n, then 



E 



4 m) (/,M)4 n) (/,M) 



0. 



It is therefore enough to establish the second equality in ( |3.2| ) since this implies ([T^), and then 
orthogonality of the terms in the series implies the first equality in ( |3.2|) . First note that, with 

Sn+l = t, 



E 



(4 n \f))'' 



t rs 



«*"" I I •••/ / f(yn+i)f(zn+i)dy n +idz n+1 G Sl fj,(yi)G Sl fj,(zi) 

JO JO Jo iR 2 ("+ 1 )< i 

n 

■ ]J [Gs i+1 -s z {yi+i - yi)G Si+1 - Si (zi+i - Zi)\yi - Zi\~ 2 dyidzidsi . (3.4) 
i=l 



Expanding the exponential in the final term of (|3.2D we have 



pt,x',y l 



exp 



o \X}-XI\ 



;ds 



1 i -pt,x',y' 

1 + ^0,x,y 



E 



i 



t rt 



n 



K dSi 



00 rt rs n 



n=l 



JO 



n=l Ul J ° J ° i=l |xj - X 2 



*2 



rpt,x',y' 
^0,x,y 



n 



i=i - x 2 . 
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Substituting this sum into the right hand side of (|3.2| ), one may match, by using the finite 
dimensional distributions of the Brownian bridge, the nth term with the expression E[lj: n \f)] 2 
in «. I 



The chaos expansion defines a linear random functional on test functions (in that there is a 
possible null set for each linear relation) . Also this linear random functional satisfies the moment 
bounds ( |1.6[ ) and ( |1.7|) . The second moment bound implies that there is a regularization (see 
I gogg ] Theorem 2.3.3), ensuring there is a random distribution ut so that 



u t ( 



(/) = E 4 n) (/> for a11 / e C c , almost surely. (3.5) 



To show that ut is actually a random measure we now construct a sequence of SPDE ap- 
proximations to ([O]). We will index our approximations by numbers e > 0. Recall that 
h(x) := \x\~ 2 = (g * g)(x), where g(x) = C7\x\~( d+2 ^ 2 . Let 

g&( x ) = (c 7 |x|-( d+2 )/ 2 ) A £~\ and h&(x) = (g& * g&) (x). 

As e I we have g^ £ \x) f g(x) and hS E \x) | /i(x). We can construct, as in (|3.1|) , a mean zero 
Gaussian field F^ e \t,x) with covariance 



F^(t,x)F^(s,y)] = 5(t - s)h ( - £ \x - y). 



We consider the approximating SPDE 

^ = I Au (e) + Ku^Fic), u (e) = /i («5) ) (3 . 6) 

with the initial condition /i^' = G$n, for some 5 = 5(e) > to be chosen later. Since the corre- 
lation is continuous in x and y, standard results give existence and uniqueness of a non-negative, 
continuous, function- valued solution u\ (t, x). Moreover we may represent the solutions in terms 
of a chaos expansion 

oo 

^ £) (/) = E4 n ' £ W 5) ) 

n=0 

where the terms I^ n ' e \f, fi^) are defined as in ( |1.16| ) and ( |1.17| ) except that fi, F are replaced 
by fj,^,F^ £ \ We now connect the approximations with the original series construction. 

Lemma 6 Suppose that /i £ TL a +- Then we may define I^ n \f, /x) and lj; n ' £ \f, H®) on the same 
probability space so that, for suitably chosen 5(e) > 0, fixed t > and f £ C c , 

u i £ \f) ~^ u t(f) in L 2 as e —> 0. 



Hence the chaos expansion (3.1) defines a random measure ut(dx), for each [i £ TC a + and t > 0. 



Proof. Let W(t, x) be a space-time white noise on [0, oo) x H d and construct both the noises F 
and F^ e ' using W as in (|3.1|). Using the convergence of both the F and F^ £ ' chaos expansions 
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and the orthogonality of multiple Wiener integrals of different orders, we find 

/•; ' 21 



(i4 B) (/) -«*(/) 

00 r 



n=0 

oo 



(3.7) 



< 2^: e (4%,M)-4 n,£) (/,M)) +2E^ (4 n %^)-4 n ' £) (/,M (5) ) 



n=0 



n=0 



We show separately that both sums on the right hand side of (^J) converge to zero as e j 0. 
We use the telescoping expansion, for n > 1, 

4 n W)-4 n,£ W) 

= K n ... / G Sl fi{yi)f(y n+ i)dy n+ i TT [G s Si (y i+ i - y;)F(c%, ds;)] 

/ .../ / G Sl (j,(y 1 )f(y n+1 )dy n+1 TT [G s s - yi)F (e \dyi,dsi 

Jo Jo JR(«+ 1 ) d ^ L 



m=l 

where Jt n ' m ' e \f) is defined to equal 



I G Sl ix(yx)f{y n+ x)dy n+ xJ{G s Si (y i+1 -yi) 

JO Jo JR(™+!) d ^ 

m— 1 n 

• [] F(dy l ,ds i ) ■ \F(dy m ,ds m ) - F^(dy m ,ds m )] ■ f[ F^(dy u d Sl ) 



i=l i=m+l 

and where a product over the empty set is defined to be 1. The isometry for the stochastic 
integral gives 

, 21 



E 



,2n 



G Sl fj,(y 1 )G Sl fj,(z )f(y n+ i)f(z n+1 )dy n+ idz n+1 

JO jR 2 ( n + 1 ) d 
n 

■ II [Gsi+i-SiiVi+i - yi)G Sl+1 - Sl {z i+ i - Zi)dyidzidsi] 



m—l 



n %* 



9~9 



i=l 



)*(g-g {£) 



{y m -z m ) \\ h^ £ \yi-Zi). 

i=m+l 



Note that < [(g - g^) * (g - < h(x) and that [{g - g&) * (g - | as 

e — > 0. Using the finiteness of E (I^ n \f,ij,)) 2 , the dominated convergence theorem implies 

. 21 



that E 



l and therefore 



lim E 

elO 



(4 n,£) (/^)- 4 n) cm 



0. 



(3.8) 
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The isometry, and h^(x) < h(x), imply that 



E 



(4%) 



Jn,s) 



(/) 



< 2E 

< 4E 



C#°(/)) s 



+ 2E 



(4 n ' £) (/) 



(3.9) 



Now, using (|3.8|) , fl3.9j), the convergence of the series Y^??=o E 



convergence theorem, the first term on the right hand side of (3.7) goes to zero as e j 0. 



and the dominated 



We now show that for fixed e > the second term on the right hand side of ( |3.7[ ) converges 
to zero as 5 1 0. Recall the initial condition was ji^ = G$n for some 5 = 5(e) > 0. But for fixed 
e the isometry shows, as in Lemma 0, that 



n=0 



5> (/ t (n,£) (/,M)-4 n,£) (/>M (£) ) 

f(y')f(x')G t (x-x')G t (y-y>) 



R 4d 



0,x,y 



exp ( jf n?h^ (Xg - X 2 )dsj dx'dy' (ji - ^) (dx) (fi - ^*>) (dy). 



When e > 0, the Brownian bridge expectation is a bounded continuous function of x,y,x',y' 
and the convergence to zero as 5 j is clear. This completes the proof of the L 2 convergence 
stated in the lemma. 



The L 2 boundedness of u t (/), for each / G C c , implies that {u^^fix} is a tight family of 
random Radon measures. The L 2 convergence of u\ (/) implies that there is a random measure 



Ut satisfying ( |3.5| ) and that u\ ^> ut in distribution as e — > 0. I 

It remains to show that {ut(dx)} is a solution of (|1.1|), and for this we must show that there 
is a continuous version of the process t — ► u% and that it satisfies the martingale problem ( |1.4| ) 
and Fix f <E C 2 . From the definition (|ITT7l ) we have, for n > 1, 

lW(y,z) = K f f G s ^(y-y')I^ 1) (y',z)F(dy',dr). 



Then using a stochastic Fubini theorem (see [Wal8£] Theorem 2.6), and the fact that Gt * f(y) 
solves the heat equation, we have, for n > 1, 



o 



1 

2 

K 
K 







lW(y,z)Af(y)v(dz)dyds 

G s . r (y - y')4 n - 1 \y',z)F(dy',dr)\ Af(y)f,(dz)dyda 
G s - r *Af(y')ds) [ 4 n - l \y',z) f ,(dz)F(dy',dr) 



kJJ (G t ^*f(y')-f(y')) J I^ 1 \y',zMdz)F(dy',dr) 
ll n \f,H)-K^ J J f(y')4 n ~ 1 \y',zMdz)F(dy',dr). 
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Rearranging the terms, we see that for each n > 1, the process 

4 n) (/) : =4 n) if^)- fl^iUf,^ 



10 '2 

is a continuous martingale. We now define 

N 



, « / / / f(y)li n ~ 1 \y,z) f i(dz)F(dy,ds) (3.10) 



A' 



UJVlt (x) = E /w^arx^). w/) = E 4 (/)• 

fc=0 fc=l 



Then, for / G and iV > 1 



/*' 1 

UN,t(f) = Kf) + / U N , S (-Af)ds + Z N ,t(f)- 
JO 2 



(3.11) 



Lemma ||] implies that E[(ujy : t(^-f) — u t(Af)) 2 ] converges monotonically to zero. Using the 
domination from Lemma [l] part i) we have 



E 


sup 


L 




t<T 





ft 1 f* 1 

/ u N:S (-Af)ds - / u s (-Af)ds 
Jo 2 JO " 



Lemma ^ also implies that ZN,t{f) converges in I? to zt(f) and by Doob's inequality 



E 



SUp 1 2jV,t(/) - z t {f)Y 
t<T 



0. 



This uniform convergence and ( 3.11 ) shows there is a continuous version of both t — > and 
t — > ut(f). Using this for a suitable countable class of C% test functions / shows that there is a 
continuous version (in the vague topology) of t — > ut- 

Now we calculate the quadratic variation zt(f), which is the L 1 limit of (zj\r r (f)) t . It is 
enough to consider the case / > 0. Using ( |3.10D we have 



{z N+1 ,(f)) t 



N+1N+1 rt 

EE* 2 " 

k=l 1=1 

rt 

.2 



(M^M jCfc' 1 ) (x, x')I^ (y, y')fi(dx')n(dy')dx dy ds 
o JR 4d \x - y\ l ' 

mm 



-u N , s (x)uN,s{y)dx dy ds 



o JR? d \x - y\ 

* f mm ( j \ u w 

/ -; r^UsldxJUsldyjds. 

lo J-R? d \x - y\ 2 

We need to justify this final convergence and we split the task into two terms 

« r mm 



(3.12) 



-UN, s {x)dx (uN, s (y)dy - u s {dy)) ds 



o Jn 2d \x - y\ 

+ [ [ f^f^ (u Ns (x)dx - u s (dx))u s (dy)ds. 
Jo Ju 2d \x - y\ z 
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We show the first term converges to zero in L , the argument for the second term is the same. 
We use the fact that \x — y\~ 2 is a convolution of c 7 \z\~^ d+2 ^ 2 with itself to see that 

;l ' " ' u NiS (x)dx (u NjS (y)dy - u s (dy)) 



r m \x-y 

C7 



mm 



C7 



R3d \ x _ z \{d+l)/2\y _ z \{d+2)/2 
UN,s(fz)(uN,s{fz) ~ U s (f z ))dz 



u NtS (x)dx(u N}S (y)dy - u s (dy))dz 



where f z (x) = f(x)\x — z\ ( rf + 2 )/ 2 . Hence, by the Cauchy-Schwartz inequality, 

rt r mm 



E 



o Jn 2d \x-y 
t 



< c 7 







E 



u NyS (x)dx(u NyS (y)dy - u s (dy)) ds 

1/2 



(3.13) 



(uN,s{fz)Y 



dzds 



E 



(.UN,s(fz) - u s(fz)Y 



dzds 



1/2 



The argument from Lemma || shows that E [(uN,s(fz)) 2 ] can be bounded uniformly in N by 



E 



[uN,s{fz)f 



< 



R 4d 



f z (x')f z {y')G s {x-x')G s {y-y')[l + 



\x — y\ a \x' — y 



)n{dx)^{dy)dx dy . 



The same bound holds for E [(u s (f z )) 2 ] ■ It is straightforward but lengthy to estimate this term. 
We show how to deal with the most singular term only. The method is to estimate the dx'dy' 
integral first using the inequalities ( |1.S| ) and ( 1.10| ). Applying Holder's inequality in the same 
way as in Lemma [l], these inequalities imply that 

mw) 



R 2d 

< 



G s (x - x')G s (y - y') - _ ^J^T^mJ^' - y ^ dx ' dy ' 

C{a)e -a\ X \-a\y\ s -a/2 ^ _ ^|-(<*+2)/2 A g -(«H-2)/4) ^ _ ^|-( d +2)/2 A ^(d+2)/^ 



< C(P,a)e~ a ^~ aM s~ 1 ~ a+( - l3/2) \x - z\'^ d+f3 ' a)/2 \y - z \-(d+/3-a)/2 . 

where we have chosen (3 £ (a, a + 2) and a so that fj, £ Tig. To apply Holder's inequality here, 
splitting the three factors f(x')f(y'), \x' — z\~^ d+2 ^ 2 \y' — z|~( d + 2 )/ 2 and \x' — y'\ a , we needed the 
bound a + + 2)/2) < d, which is implied by our assumption that a < (d— 2)/2. Substituting 
this estimate into (|3.13[ ) we find 



E 



(uN,s(fz)f 



< C(j3,a) 



dzds 



e -a\x\-a\y\ s -l+(0/2)i x _ y \ 



10 JR ad 

■\x - z\-^ +p - a ^ 2 \y - z\~ (d+(3 - a ^ 2 fi{dx)fi{dy)dzds 
= C{P,a) C [ e- a ^- a ^s- 1+( - f3 ^\x-y\- l3 fi{dx)fi{dy)dzds 

JO JR 2d 

which is finite since [i G Tip. This bound also gives the domination required to see that 

JqJE [(uN,s(fz) — Us(fz)) 2 ] dzds — > as n —* oo. This finishes the justification of the con- 
vergence in ( |3.12|) , identifying the quadratic variation (z.(f))t, and completes the construction 
of a solution {ut(dx)} to (1.1) started at fi. 
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4 Self Duality and Uniqueness 



In this section we establish the self duality of solutions in the following form: 



Proposition 1 Suppose {ut(dx)} and {vt(dx)} are solutions of with deterministic initial 

conditions uo(dx) = f(x)dx and vo(dx) = g{x)dx. Suppose also that sup x e~ a ^ f(x) < oo for 
some a and that g(x) is bounded and has compact support. Then ut(g) has the same distribution 
as v t {f). 



Remarks 

1. The duality formula is immediately clear for the solutions constructed using the chaos 
expansion in Section, || since the expression ( 1.17] ) for the nth order of the expansion is symmetric 



under the interchange of y and z. We will show in this section that the self duality relation holds 
for any solution to We then use the self duality relation to show uniqueness in law for 

solutions. 

2. Even when working with the martingale problems the self duality relation is heuristically 
clear, as can be seen by applying the technique of Markov process duality (see Ethier and Kurtz 
HEK86 ] chapter 4). Take {ut(dx)} and {vt(dx)} to be independent solutions to (IT). Suppose 



(falsely) that the solutions are function valued and have suitable behavior at infinity such that 
the integrals u s (vt- s ) and vt- s (u s ) are finite and equal by integration by parts. Take a twice 
differentiable h : [0, oo) — > R. Applying Ito's formula formally, using the martingale problem 
(|L|), leads to 

^-h(u s (v t - s )) = (l/2)h'(u s (v t - s )) (u s (Av t - s ) - v t - s (Au s )) + martingale terms. 
as 

Here we have used the cancellation of the two second derivative terms involving h" after applying 
Ito's formula for u s and for i>t_ s . Applying integration by parts the term (u s (A.vt- s ) — vt- s (Au s )) 
vanishes and this leaves only martingale terms. Taking expectations and integrating over s G 
[0, t] leads to 

E[h(u t (g))]=E[h(v t (f))} (4.1) 

which implies the self duality. To make this argument rigorous we shall argue using a smoother 
approximate duality relation. 

3. The self duality relation can be extended to hold for more general initial conditions and 
to be symmetric in the requirements on the initial conditions [i and as would be expected 
by the symmetry of the chaos expansion. One needs to define certain collision integrals (//, v) 
between measures in Ti. a+ . For example, suppose /i, v S TC a + and for simplicity suppose both 
are supported in the ball B(0,R). Define f e {x) = J 4> e {x — y)u(dy), the density of the measure 
4> £ * v. Then, if {ut(dx)} is a solutions started at fi, we claim that the variables 

ut(fe)= / / <Pe(x - y)u t {dx)v(dy) 
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are Cauchy in L 2 as e — > 0. Indeed, using the second moment formula (|1.7|), a short calculation 
leads to 



E 



(Mfe) ~ U t (fe')Y 



E 



(ut(fe-fe>)Y 



< C(t,R,fx) J J (f e (x)- f e >(x))(My)- f £ ,(y))(l + \x-y\~ a )dxdy 
= C{t,R,n)\\(j) e * v - (j) £ i * u\\ 2 a . 

Here we are extending the use of the norm \\fJ>\\ a to signed measures. Now it is not difficult 
to show that \\d> F * 

^ — ^IIq — * ^ — ^ which completes the proof of the Cauchy property. 
Denoting the L 2 limit as ut{v), and constructing vt{n) analogously, the duality relation then 
holds in this extended setting when [i, v G + , although we make no use of it in this paper. 

In the rest of this section we give the proof of Proposition |l| and deduce uniqueness in 
law and the Markov property. The proof follows from two lemmas, the first of which is an 
approximate duality relation where we smooth the measure valued solutions. 



Lemma 7 Suppose {ut(dx)} is a solution of with initial condition /x and {vt(dx)} is an 

independent solution with a compactly supported initial condition v. Suppose h : [0, oo) — ► R has 
two bounded continuous derivatives and tfi : 1Z d — ► [0, oo) is continuous with compact support. 
Fix < to < ti and a bounded o~(u s {dx) : < s < to) variable Zt . Then 



E 



Z to h 



-E 



4>(x - y)u tl (dx)u(dy) 
J t J~R 4d 



E 



Z to h 



(x - y)u to (dx)v tl - to (dy) 



(x - y)u s (dx)v tl - s (dy) <f>(xi - yi)<f>(x 2 - 2/2 



1 



1 



\X\ - x 2 \ 



\Vi ~ 2/2 1 



u s (dxi)vt 1 - s {dy 1 )u s (dx2)vt 1 - s {dy 2 )ds 



(4.2) 



Proof We first establish some integrability, sufficient to ensure that the expectations on the 
right hand side of fl4.2|) is finite. Using the independence of {ut(dx)} and {vt(dx)}, the com- 
pact support of (f) an d the bound on second moments in (|1.7| ), a lengthy but straightforward 
calculation, similar to that in Lemma [j], yields 



E 



4>{xi - yi)(f>(x 2 - y 2 ) 



+ 



u s (dxi)v t (dyi)u s (dx 2 )v t (dy 2 



/R 4d \\Xi-X2\' 1 2/1 — 2/2 1 

< C{4>^,v,T)(s-( 2 -^+l 2 +t -(2-a)+/2) for all s,t<T. 

Furthermore, using the formula for first moments (|1.6|), an easy calculation shows that 



E 



cp(x - y)u s {dx)v t {dy) 



<C((p,fi,u,T) for all < s,t < T. 



(4.3) 



(4.4) 



We now follow the standard method of duality, as explained in Ethier and Kurtz |EK86] Section 
4.4. Take / £ C 2 , apply Ito's formula using the martingale problem for ut(f) and then take 
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expectations to obtain, for s > to, 
E[Z to h(Mf))}- E[Z to h{u to {f))} 



E 



to 



Z tQ (ti(u r (f))u r (^Af) + yfc"(«r(/)) 



/( i)/( 2) ^ (dx]\u r (dx^) 
\x\ - x 2 r 



dr. 



Here Lemma [T] implies that the local martingale arising from Ito's formula is a true martingale. 
Now take : R 2rf — > R, twice continuously differentiable and with compact support. Replace 
the deterministic function f(x) by the random C 2 function, independent of {ut(dx)}, given by 
fix) = J ip(x,y)vt(dy). Fubini's theorem and the integrability in fl4.3|) and ( ]4.4[) imply that, for 
s > t , 



E 



(/. 










z to ti 


J t 






"e\ 




to 

f - 



E 



Zt h{ / / if)(x,y)u to (dx)v t (dy) 



til f ( \fj(x,y)u r (dx)v t (dy)] f f -A^ip(x, y)u r (dx)v t (dy) 



dr 



Z to h" 



tp{x,y)u r (dx)v t {dy) 
Y^-u r (dx 1 )u r (dx 2 )vt(dyi)vt(dy2] 



dr. 



\Xi — X2\ 

In a similar way, applying Ito's formula to vt(f), we obtain the decomposition 



E 



Zt h[ / / ip(x,y)u s (dx)v t (dy) 



E 



Zt h[ / / ip(x,y)u s (dx)v(dy) 



E 



z to ti 



ip(x,y)u s (dx)v r (dy)) [ f -A^ip(x,y)u s (dx)v r (dy) 



dr 



.2 r t 



+■ 



2 Jo 



E 



Z to h" 



ip(x,y)u s (dx)v r (dy) 



ip(x 1 ,y 1 )4)(x 2 ,y2) 



R 4d 



\yi - V2\ 



Now defining 



F(s,t) = E 



Z to h 



u s (dxi)u s (dx 2 )vr (dyi )v r (dy2 
ip(x,y)u s (dx)v t (dy) 



dr. 



the last two decompositions show that s — > F(s, t) and t — > F(s, t) are both absolutely continuous 
and gives expressions for their derivatives d\F(s,t) and d2F(s,t). Then applying Lemma 4.4.10 
from [EK86| we obtain 



E 



Zt ° h \J J $( x >y)^i( dx ) v ( d v) 

-- F(t 1 ,0)-F(t ,t 1 -t ) 



E 



Zt h[ / / f(x,y)u to (dx)v tl -t (dy) 



d\F(s, t\ — s) — d2F(s, t\ — s)ds 

ti— to 
h 

E 

to 

2 rti 

i/}(x,y)u s {dx)v tl - s (dy) ) rp(xi,yi)ip(x 2 ,y2) 



Z to ti (J J i>{x,y)u s {dx)v tl - s {dy)\ J J \{& {x) ~ A^(x,y)u s (dx)v t ^ s (dy) 



ds 



K 



2 Jto JR id 



E 



Zt^ti' 



1 



1 



\x\ -x 2 \ 2 I2/1 — 2/2 1 2 



u s (dxi)v tl -s{dyi)u s (dx 2 )v tl - s {dy 2 ) 



ds. 



(4.5) 
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Now suppose that (j) : R d — > [0, oo) is smooth and has compact support. Choose a series 
of smooth, compactly support functions ip n (x,y) satisfying 0<^ n |lasn->oo and with 
d x ip n , dyipni d xx ip n , dyyipn converging uniformly to zero. Apply (|4.5| ) to the function ifj(x,y) = 
il) n (x,y)<t>{x - y). Using (A^ - A(w))^(a: - y) = we may, using the integrability in ( [4.3| ) and 
(4.4), pass to the limit in (|4.5| ) to yield ( |4.2| ). Finally we obtain the result for general continuous 
4> by taking smooth approximations. I 



Now we take 4>(x) a smooth, non- negative function on R d , supported on the unit ball 
{x 6 : \x\ < 1} and satisfying J- Rd (j)(x)dx = 1. Define an approximate identity by <fi e (x) = 
e~ d 4>{x / s). We may and shall suppose that < 4>(x) < 2G £ (x) and hence that 4> £ < G £ 2. We 
shall use this test function in Lemma ^) and the following lemma controls the right hand side 
of (P). 



Lemma 8 Suppose {ut(dx)} and {vt(dx)} are independent solutions of (1.1), with initial con- 
ditions fjL, v , where v compactly supported. Then 



E 



lo Jn 4d 
converges to zero as e — > 0. 



1 



1 



\x\ - x 2 \ 



m - V2\ 



u s (dxi)vt-s(dyi)u s (dx2)vt- s (dy 2 )ds 



Proof This lemma is a straightforward but lengthy consequence of the second moment bounds 
(|1.7|). Since it is this proof that requires the strict inequality k < (d — 2)/2 and also the 
requirement that /i, u € Hp for some (3 > a, we give some of the details. 



The second moment bounds show that show that the expectation in the statement of the 
lemma is bounded by 

C j j ^{x'y - y'i)</> 6 (x 2 - y' 2 )G s (xi - x' 1 )G s (x 2 - x 2 )G t - s (yi - y' 1 )G t - s (y 2 - y' 2 ) 

Jo JR M 

x , (t-s) a 



1 + 



\X\ — X2| a |x' 1 — X 2 \ a 



1 + 



vi -y2\ a \y[ -y' 2 \ a 



I 12 



\y\ - y'2 



' 12 



^{dxi)^{dx 2 )u{dyi)v{dy 2 )dx' l dx 2 dy'idy' 2 ds. (4.6) 



The idea is to bound first the dx^dx^dy'idy^ integral. We can split the dx^dx^dy^dy^ integral 
into four terms by expanding the brackets 



1 + 



\%1 — X2| a |x' 1 — X 2 \ a 



1 + 



(t - s y 



vi -y2\ a \y[ -y' 2 \ a 



We shall show only how to treat the worst of these terms, namely 

, (t>e{x'i - y'x)<f> s (x 2 - y' 2 )G s (x 1 - x' 1 )G s (x 2 - x' 2 )G t - s (yi - y' 1 )G t ^ s (y2 - y'2) 
s a (t-s) a 



R 4d 



(4.7) 



\xx - x 2 | a |x' 1 - x 2 \ a \yi - y 2 \ a \y[ - y 2 



1 \a 



1 



1 



Wx~ x 2? \y'i-y'2? 



dx 1 dx 2 dy l dy 2 ds. 



This is the term that requires the restriction on n. The other three terms are similar but easier. 
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We split the integral (|4.7| ) into two regions. First we consider x'i,y'i,x' 2 ,y 2 i n the set 

A £ = {\x'i — x' 2 \ > £ 7 , \y'i — y' 2 \ > £ 7 }, where 7 G (0, 1) will be chosen later in the proof. On this 
set, since we may also suppose — y'-y\ < e and |x' 2 — 7/2 1 — e by the support of <j> s , we have, 
arguing using the mean value theorem, 



1 



I / / 1 9 

I t C 1 n 



/ 12 



< C(7)e 1 ^ 37 for all e < 1/3. 



a J/ 2l 1 2/i 

This bound means the integral ([Ol ), over the set A e , can be bounded by 

s a {t - s) 



C(j)e 



1— 37— 2«7 



\xi - x 2 \ a \yi - y 2 \ 



-G t+£ 2{xi - yi)G t+£ 2{x 2 - y 2 ). 



We shall choose 7 > so that 1 — 37 — 2cry > 0. It is easy to show that this bound substituted 
into Q4.6| ) will vanish as e j 0. To estimate the integral fl4.7| ) over the complimentary set A £ we 
simply bound 



1 



I ' ' |2 



/ 12 



< 



1 



1 12 



+ 



1 



7 12 



and it becomes 

s a (t- s) c 



G s (xi - x[)G s (x2 - x' 2 )Gt-s{yi - y[)G t -s{y2 - y 2 ) 



\xi - x 2 \ a \yi - y 2 \ a ja% 

■4>ei x 'i ~ Vi)4>e{x' 2 - y' 2 )\x' x - x' 2 \~ a \y' x - y' 2 \~ a (\x[ - x' 2 \~ 2 + \y[ - y' 2 \~ 2 ^ dx^dx^dy^dy^ 

^ ^ I G s {xi - x' l )G s (x 2 - x' 2 )G t -s(yi - y' 1 )G t - s (y 2 - y' 2 ) 

xi - x 2 \ a \yi - y 2 \ a Ja% 

■Mx[ - y[)M4 - y' 2 ) (\x[ - x' 2 |-( 2+2a ) + \y' x - y' 2 \^ 2+2 ^) dx' x dx' 2 dy' x dy' 2 . 



< C 



We show how to treat just the integral with the term \x x — x' 2 \~( 2+2a \ the term \y[ — y' 2 \~^ 2+2 °^ 
being entirely similar. Note that the restriction k < {d — 2)/2 is simply to ensure that 2 + 2a < d 
and so the pole |,z| _ ( 2+2a ) is integrable on H d . We may choose 5 £ (2 + 2a, (2 + a + j3) A d). 
Then, using the bound cj> e < 2G £ 2, we can do the dy[dy 2 integrals to see that 



G s {xi - x[)G s (x 2 - x' 2 )G t -s{yi - y'i)G t -s{y 2 - y 2 ) 

A* 

■<j>e(x x - yi)<j>e{x 2 -y' 2 )\x'i- x' 2 \~ i2+2a) dx[dx' 2 dy' x dy 2 

< / G s (x x - x[)G s (x 2 - x' 2 ) 

J{\x[-x' 2 \<e} 

■G t - s +e 2 (yi ~ Xi)Gt- s + e i(m - x' 2 )\x[ - x' 2 \~( 2+2a > dx' x dx' 2 

< Ce^ 5 - 2 - 2 ^ [ G.(xi - x[)G s (x 2 - x' 2 ) 

■G t _ s+E 2(yx - x' 1 )G t _ s+£ 2(y 2 - x' 2 )\x' x - x' 2 \~ s dx' x dx' 2 . 
We now split into two cases: s < t/2 and s > t/2. When s < t/2 we have the bound 

G t _ s+£ a(yi - x[)G t _ s+£ 2(y 2 - x 2 ) < C(a,u,t) ex.p(-a\x[\ - a\x' 2 \), for y x ,y 2 Gsupp(» 



2S 



So, when s < t/2, 

G s {xi - x' 1 )G s (x 2 - x' 2 )G t -s+e{yi - x' 1 )G t - s+ e{y2 - _ x 2 \~ 5 dx^dx'z 

< C(a,u,t) / G s {x\ — x' 1 )G s (x2 — x'2) expf— a\x'i\ — a\x' 2 \)\x'i — x' 2 \~ s dx'idx 2 

< C(a, u, t) exp(— a\x\\ — a|x2|) [\x\ — x 2 \~ 5 A s~ s ^ 2 ^j 

< C{a,u,t)^ V (-a\x l \-a\x 2 \)\x l -X2\- t3+a s- {& -f i+a)l2 , 

using the tricks from Lemma [l] for this last inequality. Combining all these bounds one has, 
when substituting the integral ( [4 .7] ) over the region A c e into (4.6), and considering only the time 
interval [0, t/2], the estimate 

C(a,u,t)e (5 - 2 - 2a) ^ / / {-. m . ]e- a ^- a \ X2 ^(dx l )fi(dx2)u(dy 1 )u(dy 2 )ds. 



,0 JR 4d \ \xi - x 2 | /3 |yi - y 2 \ a ; 

Choosing a so that e~ a ^ fi(dx) G TC®, the integral is finite and so this expression vanishes as 
e I 0. The integral over [t/2,t] is treated in a fairly similar way using the assumption that 
v G Tip. I 

To deduce Proposition [l| fr om Lemmas [7] and || is easy. By a simple approximation argu- 
ment it is enough to prove ( |4.1| ) for h with two bounded continuous derivatives. We apply the 
approximate duality relation ( |4.2| ), using = to < t\ = t and Zt a = 1, to the function <f> e . Then 
take e — > and use the control on the error term in Lemma § to obtain the result. 

We show two consequences of the duality relation and its proof. 



Corollary 1 Solutions to are unique in law and we let denote the law of solutions 

started at fj, £ Tt a + ■ 

Proof First suppose that {ut(dx)} and {vt(dx)} are two solutions with the same deterministic 
initial condition /i. Construct a third solution {wt(dx)}, independent of {ut(dx)} and {vt(dx)} 
and with initial condition wo(dx) = f(x)dx for some non-negative, continuous, compactly sup- 
ported function /. Then apply the approximate duality relation ([4.2|), with = to < t\ = t and 
Z to = 1, to the pair {ut(dx)} and {wt(dx)} and to the pair {vt(dx)} and {wt(dx)}, using the 
function cj) £ . Subtracting the two approximate duality relations we see that 



E 



M / / &e( x ~ V) u t{dx)f{y)dy 



E 



h [ 4>e{x - y)v t {dx)f{y)dy 



equals the sum of two error terms, both of which converge to zero as e — > by Lemma ||. Hence 
E[h(ut (/))] = E[h(vt (/))] for all such / and for all suitable h. Choosing h(z) = exp(— Xz) we 
obtain equality of the Laplace functionals of ut(dx) and vt{dx) and hence equality of the one 
dimensional distributions. 

Now we use an induction argument to show that the finite dimensional distributions agree. 
Suppose the n-dimensional distributions have been shown to agree. Choose < s\ < s 2 . . . < 
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s n+ \ and set t\ = s n+ i,to = s n . Then apply the approximate duality relation ( |4~2^ ) to the 
pair {ut(dx)} and {wt(dx)} with Z to = 117=1 ex P( — u s t (fa)) for compactly supported fa > 0. 
Also apply the approximate duality relation (|4.2|) to the pair {u t (cZx)} and {w t (dx)} with 
Z io = nr=i ex P(~ v Si (fa))- Subtracting the two approximate duality relations, use the equal- 
ity of the n-dimensional distributions and let e j to obtain equality of the n + 1-dimensional 
distributions, completing the induction. Since the paths have continuous paths the finite dimen- 
sional distributions determine the law. 



For general initial conditions uq we let P^ be a regular conditional probability given that 
uo = fx. It is not difficult to check that for almost all /i (with respect to the law of uo) the process 
{ut(dx)} is a solution to (1.1) started at \i under P^, (The moment conditions carry over under 
the regular conditional probability and these allow one to reduce to a countable family of test 
functions in the martingale problem). By the argument above the law of {ut(dx)} under the 
conditional probability P^ is uniquely determined (for almost all fi). This in turn determines 
the law of {ut(dx)}. I 



Corollary 2 For any bounded Borel measurable H : C([0, oo), M) — ► R the map \x — ► Q^[H], 
the integral of H with respect to Q^, is measurable from 7i. a + to R. 

The set of laws {Q^ : \i E T~C a +} forms a Markov family, in that for any solution {ut(dx)} 
to for any bounded measurable H : C([0, oo), M) — > R, and for any t > 

E [H(ut+-)\J : t\ = Qu t [H] i almost surely. 



Proof. We use the methods of Theorem 4.4.2 of Ethier and Kurtz [ EK86 ]. We were unable to 
directly apply these results, but with a little adjustment the methods apply to our case and we 
point out the key changes needed. 



We only allow initial conditions in the strict subset TC Q + of all Radon measures, and do not 
yet know that the process takes values in this subset. But by restricting to the ordinary Markov 
property it is enough to know that P(ut(dx) E T~l a +) = 1 for each fixed t, and this follows from 
Lemma |l] part ii). 



The measurability of fj, — * Qu\H\ can often be established for martingale problems by 



establishing it as the inverse of a suitable Borel bijection (see [EK86] Theorem 4.4.6). We do 
not use this method as TC a + is not complete under the vague topology. However the measurability 
can be established directly as follows. It is enough, by a monotone class argument, to consider 
H of the form H(oS) = n?=i hi( u ti(fa)) fo r bounded continuous functions hi, for fa E C c , for 



o < h < t 2 . . 

from Section 



t n and for n > 1. But for such H we can write, using the construction of solutions 



IIM£4 n W) 



U=l 



\n=0 



lim E 

N~>oo 



n ^ 



N 



\n=0 



For each N < oo the integrands J2n=o (f^ I 1 ) are > °y * ne definition of the maps I^ n \f,n), 
continuous in \x. So Q^[H] is the limit of continuous maps on 7i a +- 



30 



We can now follow the method of in Theorem 4.4.2 part c) in Ethier and Kurtz | EK86 | in 
the proof of the Markov property. The only important change in the argument from Ethier and 
Kurtz is that we have uniqueness in law for solutions to (|1.1|), and this requires the moment 
bounds ( |1.6| ) and ( |1.7| ) to hold as well the martingale problem ( |1.4| ) and (1.5). The key point is 
to show that, for any t > 0, the process {ti t+ .(dx)} satisfies these moment bounds. For this it is 
enough to show for all / : R rf — > [0, oo) and < s < t 



E[ut{f)\a{u r {dx):r<s)] 



G t - S (x - x')f(x')u s (dx) 



and there exists C, depending only on the dimension d and k, so that 

2 



E 



f(x)u t (dx) 



a(u r (dx) : r < s) 



C 



R 4d 



G t - S (x - x')G t - s {y - y')f{x')f{y') 1 + 



\x — y\ a \x' 



u s (dx)u s (dy)dx' dy' . 



By uniqueness in law it is enough to prove these bounds for the solutions constructed via chaos 
expansions in Section ||. It is also enough to prove these bounds for / 6 C c . The first moment 



follows from the fact that E[l\ ■ (/, ju)|.Fs] = Is' 1 ' (Gt- S f, /•*) and the convergence of the series 
(|3.2| ). For the second moment bound we use the approximations introduced in Section ^, for 
which we know u^- e \f) — ► u t (f) in I? . Fix < s\ < . . . < s n < s, fx, . . . , f n 6 C c and a bounded 
continuous function h : R" — > R. Then, using the Markov property of the approximations u 



(0 
t i 



E 



( Ut (f)) 2 h(u Sl (f 1 ),...,U Sn (fn)) 



limE 



< ClimE 

elO 



(u ( f\f)) 2 h( Usi (f 1 ),...,U Sn (fn)) 

G t - S (x - x')G t -s{y - y')f{x')f{y) 

u^{dx)u^{dy)dx'dy' h{u si {f x ),. 



R 4d 



1 + 



(t - s y 



x — y\ a \x' 



,U Sn (fn)) 



CE 



R 4d 



G t -s{x - x')G t - s (y - y')f(x')f(y') 



(4. 



1 + 



it - s y 



x — y\ a \x' 



V 



u s (dx)u s (dy)dx' dy' h{u Sl (f\), . . .,u Sn (f n )) 



The last equality follows by the convergence vS £ '{f) — > Ut(f) for compactly supported / and 
an approximation argument using the uniform second moment bounds on and u s . The 
inequality (4.8) implies the desired second moment bound and completes the proof. I 



5 Death of solutions 



We adapt a method from the particle systems literature to study questions of extinction. Liggett 
and Spitzer used this technique, described in Chapter IX, Section 4 of Lig_8J], to study analo- 
gous questions for linear particle systems. The corresponding result for certain linear particle 
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systems, indexed on Z and with noise that is white in space, is that death of solutions occurs 
in dimensions d = 1, 2 for all k, and in dimensions d > 3 for sufficiently large k. The long range 
correlations of our noise lead to different behavior, an increased chance of death, and death 
occurs for all the values of d > 3 and k that we are considering. However our basic estimate in 
the proof of Proposition [2] below leaves open the possibility that the death is extremely slow. 

We start by considering initial conditions with finite total mass. To study the evolution of 
the total mass we want to use the test function / = 1 in the martingale problem. The next 
lemma shows this is possible by approximating / by suitable compact support test functions. 



Lemma 9 Suppose that {ut(dx)} is a solution to started at fi & Tt^- Then the total mass 
{vk(l) : t > 0} is a continuous martingale with 

{u{1))t =!J J \ x -vv ds - 



Proof. We first check that the assumptions on the initial condition imply that E[ut(l) 2 ] < oo. 
The bound on second moments (|1.7|) implies that 



E 



u t (l) 2 < C G t (x-x')G t (y-y')(l + 



< C 



R2d y \x — y\ 



t a 



\x — ?/|"|ar — y 
fi(dx) fj,(dy) using (1.9) 



, _ na dx dy n{dx){i{dy) 



(5.1) 



We may find /„ e C c 2 (R d ) so that < /„ | 1 and \\&f n \\oo i as n — > oo. Applying Doob's 
inequality we have, for any T > 0, 



E 



sup|z t (/ n ) - Zt(fm)\' 
t<T 



< CE 
= CE 

< C\E 



ZT(fn) ~ Z T {f m )\ 2 

1 f 1 

UT{fn - fm) ~ Kfn ~ fm) ~ ~ / U s (Af n - Af m )ds 

2 JO 



U T (fn ~ fmf\ + Kfn ~ fm? + ^(WJoo + II &fm\\oo)E 



T 



u s (l)ds 



This expression is seen to converge to zero as n, m — * oo by using dominated convergence and the 
bound in ( |5.1| ). From this we can deduce that, along a subsequence, zt(f n ) converges uniformly 
on compacts to a continuous martingale. Also 



Since 









2" 




E 


sup 


[ u s (Af n )ds 








t<T 


Jo 




Jo 



Zt(fn) + I U s (-Af n ) = U t (fn) - fi(fn) -> u t (l) - 
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we can conclude that «t(l) is a continuous martingale. Moreover we claim that 

u s (dx)u s (dy) 









E 


sup 


a 




t<T 





E 



\x-y 

1 - fn{x)fm{y) 



2 ds - (z(f n )) t 

u s (dx)u s (dy)ds 



/o ./ ./ \x-y 
This follows by dominated convergence and the bound 

u s (dx) u s (dy) 



as n — > oo. 



E 



-ds 



\x -y 



lim E{z{f n )) T 



lim E 

n— >oo 



(5.2) 



UtUu) ~ fi(fn) ~ u s (Af n )ds 

= E[u T (l) 2 } -/x(l) 2 < oo. 

Using fl5.2j ) it is now straightforward to identify the quadratic variation of as in the state- 
ment of the lemma. I 

Proposition 2 Suppose that {ut(dx)} is a solution to ( [Op started at (J, £ Ti a + n 7i^. Then 

lim Ut(l) = almost surely. 



Proof. The previous lemma shows that the process ut(l) is a non-negative martingale and 
hence converges almost surely. We will show that 



lim E 



Ut 



(1) 1/2 







(5.3) 



which then implies that the limit of «t(l) must be zero. We consider first the case that fi is 
compactly supported inside the ball B(0,K). We let Ct = B(0,Rt) be the closed ball with 
radius 

^ = K+[c 8 (tVl)loglog(tV4)] 1/2 

where c$ is a fixed constant satisfying eg > 4. We write for the complement of this ball. Let 
To be the first time t > that Ut(l) = 0. (In a later section we shall show that P(tq = oo) = 1 
whenever /i(l) > but we do not need to assume this here.) Using Ito's formula, and labeling 
any local martingale terms by dM, we find that for t < tq, 



Ait(l) 1 / 2 = dM t -^ t (l)V2 



u t {dx)u t (dy) 



dt 



< dM t 
= dM t 

< dM t 



2,2R 2 t 



32R 2 



32R? 



«t(l) 



n t (l) 2 |x - y\ 2 
i 2/ f u t {dx)ut{dy) 



(5.4) 



1/2 



C t JCt 

1 - 



utiV) 2 



dt 



dt 



ut 



ut(dx 
c c t «t(l) 
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< dM f 



dM t 



32i? t 2 



32R? 



u t (dx) 



Ml) 1/2 1-2 , M1 



1/2N 



dt 



The local martingale term in (5.4) is given by dMt = (l/2)u$ (l) _1//2 <itij(l) and is reduced by 
the stopping times T\i n = inf{t : ut(l) < l/n}. So applying (|5.4| ) at the time t/\Ti/ n and taking 
expectations we obtain 



E[ut An/n (l) 1 ^] < M (i)^ _ £ 



1/2 



*AT 1/n K 2 







32i?2 



■u s (l) 1/2 (fs 





7 




Jo 



u s {C c s ) 1/2 ds 



Letting n — > oo, using monotone convergence and the moments established in ( |5.lD , we obtain 
the same inequality with T\i n replaced by To- Since the paths of a non- negative local martingale 
must remain at zero after hitting zero we may further replace t A tq by t in the inequality. 
Defining rj t = E ^(l) 1 / 2 



we therefore have 



Vt < Vo 



* k 2 



o 32i? 2 



r] s ds + 



* K 2 



E 



u s (C 



CNl/2 



ds 



(5.5) 



The aim is to estimate the expectation in this inequality and to show that it implies that Vt — ► 0. 
Let 

H(s ' t) = exp {-f s im dr )- 

It follows from the definition of Rt that ^^zdr = oo for any s > and so, for any s > 0, 

lim E(s,f) = 0. (5.6) 

t— >oo 

Applying Gronwall's inequality to (pTq), we obtain 

Vt < mt)Vo + J^(s,t)-^E [usiC^ 2 ] ds. (5.7) 

If we show that / °° j^jE u s (Cf) 1 / 2 ds < oo it then follows that Vt j as t — ► oo (use 
< H(s,t) < 1, (|5.6|) and dominated convergence). Using the Cauchy- Schwartz inequality and 
the formula for first moments, we obtain 



(E [ut(Ct) 1/2 ]) 2 < E[u t (Ct)) 



< 



C? J~R d 



G t (x - y)n(dx)dy 



< Cn(l) / (27ri)- d/2 exp(-r 2 /2t)r d - 1 dr 

jRt-K 

= C/x(l) / exp(-s 2 /2)s d - 1 (is 

J(Rt-K)/Vi 

< C/x(l) exp(- (i? t - K) 2 /2t) (l + (i2( - K/V~t) d - 1 

< C/i(l) [log(t V 4)]~ cs/2 [loglog(t V 4)] (d ~ 1)/2 . 



(5.8) 
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Here we have used the following standard inequality: by the change of variables y = x + z we 
find 

/*oo /*oo 

/ ^ V {-y 2 /2)y d - 1 dy < Cexp(-a; 2 /2) / exp(-z 2 /2)(^ _1 + z d ~ l )dz 

Jx JO 

< C(l + exp(-x 2 /2), when x > 0. 
Finally we use (|5.8|) to derive the following: 



oo K 2 



IQR 2 



E[u s (C c s )] 1/2 ds<C(K,fi) 1 + 



00 (log log s)^ 5 )/ 4 



s(log s 



,cs/4 



cis < 00. 



This completes the proof in the case fi is compactly supported. In the general case we fix e > 
and split the initial condition so that \x = \jP^ + \j^> where < e and H<® is compactly 

supported. By uniqueness in law we may consider any solution with initial condition [i and we 
choose to construct one as follows: let be the strong solutions, as constructed in Section 

3, with respect to the same noise and with initial conditions /xW,//( 2 ) and set u = + u^ 2 \ 
It is easy to check that u is a solution starting at /i, which is a statement of the linearity of the 
equation. Using Cauchy-Schwartz and the formula for first moments ( |l.6|) we have 



E 



Mi) 



1/2 



= E 
< E 



n( 1) (l) + n| 2) (l)) 



1/2' 



II 



< E 
= e^ + E 



1} i/2 
1/2 



4 2) (1) 1/2 



,( 2 )/ 



«r(i) 1/a 



H 2) (i) i/2 ; . 

Thus (|5.3[) follows from the compactly supported case and the proposition is proved. 



Proof of Theorem Firstly the case of an initial condition with finite total mass. If 
P(u e H° a+ ) = 1 then, 

P(u t (l) - 0) = / Q„(U t (l) -> 0)P(u G dM) = 1- 

Secondly the case where of an initial condition that has locally bounded intensity. For such uq 
we have, using the first moment formula, 



E[ Ul (cf>)}=E 



G\(x — z)4>(x)uo(dz)dx I < C / 4>(x)dx 



(5.9) 



for some constant C < 00. That is -E[ui(<ix)] < CL(dx) where we write L{dx) for Lebesgue 
measure. Fix a bounded set A. By the linearity of the equation the map [i — » Q^(Ut(A) A 1) is 
increasing in ^u. Moreover it is concave in /i. Indeed if u^(dx) and u\{dx) are solutions started 
from /x and z^, with respect to the same noise, then, by linearity and the concavity of f(z) = z Al, 

Qe^(i-e)u(Ut(A) A 1) = E [(9u?(A) + (1 - 0)<(A)) A 1] 

> £7[(9«(A)A1) + (1-6»)«(A)A1)] 
= OQ^UtiA) A 1) + (1 - 9)Q v (U t (A) A 1). 
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Then 

E[ut+i(A) A 1] = / Qfj,(Ut(A) A l)P{u\ G (f/i) (by the Markov property) 

< Q^Md*')]^- 4 ) A 1) (by Jensen's inequality) 

< QcL(dx)(U t (A) Al) (using Q) 
= Qj A (Ct/i(l)Al) (by self duality) 

which converges to zero by Proposition |2[ This completes the proof of Theorem ^. I 



6 Support Properties 

In this section we establish the various properties listed in Theorem ^. 



6.1 Dimension of Support 



We can apply Frostman's Lemma (see | F85 | Corollary 6.6) to obtain a lower bound on the 
Hausdorff dimension of supporting sets for solutions ut(dx). Indeed Lemma [l] part ii) and 
Frostman's Lemma imply that any non-empty Borel supporting set for the measure Ut(dx), at 
a fixed t > 0, must, almost surely, have dimension at least d — a. We prove in Subsection 
that if fi ^ then ut 7^ almost surely. This establishes the fixed t result in Theorem |3] i) . We 
now show a weaker lower bound that holds at all times. 



Proposition 3 Suppose that {ut(dx)} is a solution to (1.1). 

i.) If P{uq G TL(cL-2~a)-) = 1 then P{u t G Ttu-2-a)- f or a M t > 0) = 1. Indeed, for some a, 

\p < oo J =1 for all t > and (5 < d — 2 — a. 



P ( There exists a so that sup||^e a ' s '| 



ii.) For any initial condition we have P(ut G Tlu-2-a)- f or a ^ t > 0) = 1. 



Remarks 

1. Since T-Lu-2-a)- ^ (which requires k < (d — 2)/2)) we also have, for any initial 
condition, P(ut G TC a + for all t > 0) = 1. 

2. Using Frostman's Lemma, part ii) of this proposition implies that, at all times t > 0, a 
Borel set At that supports ut{dx) must have Hausdorff dimension at least d — a — 2. 

3. The idea behind the proof of Proposition |3| is to show, for suitable values of p, that 
the process of = J f ut(dx)ut(dy)/\x — y\ p is a non-negative supermartingale. Applying Ito's 
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formula formally, ignoring the singularity in \x — y\ p , and writing dM for any local martingale 
terms, we find 

dS { t p) = dM + J J u t (dx)u t (dy) Qa (\x - y\' p ) + k 2 \x - y\- {p+2) j dt 

= dM + (p 2 - (d - 2)p + k 2 ) /" /" u t (da>t(dy)|x - y|^ +2) <ii 

where A is the Laplacian on R 2d , acting on both variables x and y. The solution to the inequality 
p 2 — [d — 2)p + k? < gives the condition a < p < d — 2 — a. The rigorous calculation below 
does not quite apply to the boundary value of p = d — 2 — a. 

First we prove a lemma extending the martingale problem to test functions on H 2d . 



Lemma 10 Suppose that {ut(dx)} is a solution to with initial condition p. Then for twice 
differentiable function f : R 2rf — > R with compact support 

M t (f)= J J f(x,y)u t {dx)u t (dy)- fj J (^Af(x,y) + « 2 /^) u s {dx)u s {dy)ds (6.1) 

defines a continuous local martingale. 



Proof. For / of product form, that is f(x,y) = J2k=l c Pk(x)'4'k(y) where 4>k^k £ C 2 , this 
claim is a consequence of the martingale problem ( |1.4|) and ( |1.5[) together with integration by 
parts. Now we claim that we can choose f n (x, y) of product form, and with a common compact 
support, so that f n and Af n converge uniformly to / and A/. One way to see this is consider 
the one point compactification E of the open box {(x,y) : \x\, \y\ < N} and to let (X t ,Y t ) be 
independent d-dimensional Brownian motions absorbed on hitting the boundary point of E. 
Then consider the algebra A generated by the constant functions and the product functions 
(f>(x)ip(y), where (j>,ij) are compactly supported in {x : \x\ < N}. The Stone- Weierstrass theorem 
shows that this algebra is dense in the space of continuous functions on E and the transition 



semigroup {T t } of (X t , Y t ) maps A to itself. A lemma of Watanabe (see [ EK86 ] Proposition 3.3) 



now implies that A is a core for the generator of (Xt,Yt) and this implies the above claim. 

The continuity of t — > ut, and the calculation in Lemma |l| part ii), imply that Mt(f n ) 
converges to Mt(f) uniformly on compacts, in probability. So the limit Mt{f) has continuous 
paths. Also, if f n and / are supported in the compact set A, the stopping times 

T k = M{t : u t (A) +f I I U ^ )U f y) ds > k} 
Jo Jaja \x - y\ z 

satisfy Tj, f oo and reduce all the local martingales Mj(/ ra ) to bounded martingales. We may 
then pass to the limit as n — > oo to see that Mt(f) is a local martingale reduced by {Tk}. I 

Proof of Proposition ^|. For part i) we may, by conditioning on the initial condition, suppose 
that uq = p £ H(d_2- a )-- We may then choose a so that p £ TL a p for all (3 < d — 2 — a. 

We shall approximate \x — y\~ p by a sequence of compactly supported functions as follows. 
Choose <p n £ C 2 satisfying and 4> n (x) = 1 for \x\ < n and with <j> n , d Xi 4> n , d XiX ,<p n uniformly 
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bounded over x and n. Define, for e > and p € [(d — 2)/2, ci — 2 — a], 

fn*(x, y) = (l + (e + \x- yfy^) e-^+N^-^+I^^^^^O/). 
A calculation shows that 

(e + |s - I/I 2 )""/ 2 



2 |x — y 

- ( 

< I e + \x — y\ 1 £ k |x — y| 



.: + |s - y| 2 ) (P+4)/2 ((p 2 - (d - 2)p + k 2 )|x - y\ 2 + (2k 2 - pd)e + e 2 K 2 \x - y\~ 2 



< C(p)\x-y\- {p+2 l 

The penultimate inequality follows from the restriction on the value of p and k and the last 
inequality follows by considering separately the cases e < \x — y\ 2 and e > \x — y\ 2 . A simple 
calculation also shows that 

A (V«(i+M 2 ) 1/2 -a(i+M 2 ) 1/2 )| < C 7( a ) e -«(i+NI 2 ) 1/2 -«(i+|y| 2 ) 1/2 < C{a)e- a \ x \- a \ y \. 
Now a lengthy calculation, using the above two bounds as key steps, shows that 



\x - y\ 



fnA x >y) 



< C(a,p) (l + \x-y\ 



(P+2)\ p -a\x\-a\y\ 



Note the bound is uniform over n and e. Using the test function fn tE (x, y) in Lemma 1C we have 
that 

M t (/ n>e ) = y J f n , e (x,y)ut(dx)u t (dy) - £(e,t) (6.2) 
is a continuous local martingale and 



£(e,t)<C(a,p) f f f (l + \x-y\ 

JO JR d JR d v 

Now we apply Doob's inequality in the following form 



-(P+2H „-a\x\-a\ : 



'u s (dx)u s (dy)ds. (6-3) 



Lemma 11 Suppose {At}, {Mt}, {Dt} are continuous processes satisfying < At = Mt + Dt 
and where Mt is a continuous local martingale with Mq bounded. Then for A > 

P (sup A, > 2A ] < ~ [E[M ] + 3E[sup\D s \]) . 

\ s<t I A \ s <t I 



Proof If {Tk} reduce the local martingale Mt then by Doob's inequality for positive submartin- 
gales 

P( sup \M S \ > A) < ~E[\M tATk \] 

s<tAT k A 

< ^(E[A tATK ]+E[\D tATK \]) 

< j(E[M ] + 2E[\D tA T K \}) 

< Ue[M q ] + 2E[sw\Ds\]\. 
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Let k — > oo and combine with the bound P(sup s<i |D s | > A) < E[sup s<t \D s \]/ X to complete 
the lemma. I 



We apply this lemma to the decomposition (|6,2[) together with the bound Q6,3|) to obtain 




(l + \x — y\ 



lim P sup 

e-»0,n-»oo \ s<< 



-a|x|— a|; 



'ut(dx)ut(dy) > 2 A 



f n , £ (x,y)u t (dx)u t (dy) > 2A 



< 



( U ''' p) -||Mdx)exp(-a|x|)||2 



A 



+ 



/oV/M^ 1 



■(p+2) 



-a|x|— a|; 



'u s (dx)u s (dy)ds 



(6.4) 



A little effort, as in Lemma [T] part ii) and using the fact that \i G "H^ for all (3 < d — 2 — a, 
shows that the expectation on the right hand side of ( |6.4| ) is finite. One needs, however, the 
strict inequality p<d — 2 — a so that the worst pole is \x' — y'\~(P +2+a \ which is therefore 
still integrable ensuring the bound ( |1.12| ) applies. The bound in ( |6.4[ ) implies part i) of the 
Proposition. 

For part ii) we may suppose, by conditioning on the initial condition, that uq = \i £ Ti-a+- 
But then Lemma || part ii) implies, for fixed to > 0, that ut (dx) £ Hd-2-a almost surely. The 
Markov property of solutions and part i) then imply that the desired conclusion holds for t > to- 
Letting to j completes this proof. I 



Corollary 3 The family {Q^ : fj, G T~C(d-2-a)-} ^ s a strong Markov family. 



Proof. Let {ut} be a solution defined on (O, J 7 , F t , P) and satisfying P(uq £ Tifd-2-a)-) = 1- 
Let r < oo be a Tt stopping time and r(n) < oo be discrete stopping times satisfying t„ j r. 
Fix < ti < . . . < t n and /i,...,/ n eC c and set H(u) = exp(i(u tl (/i) + . . . + u t „(f n )))- Fix a 
set A € J- T . Then the ordinary Markov property implies that 



E 



%n)f)l(A) =E Q Ut( JH}1(A) 



(6.5) 



If we can pass to the limit as n — > oo to replace r(n) by r, then this identity implies the result. 
By the continuity of paths the left hand side of ( |6.5| ) converges as desired. We claim that 

if \x n — > fi vaguely and sup n \\/j, n (dx)e~ a ^ \\ a < oo then Q^H] — * Q^[H]. 

Assuming this claim, Proposition ^ part i) allows us to pass to the limit on the right hand side of 
(|6.5|). To prove the claim we let ut(dx) be the solution starting at fi constructed using the chaos 
expansion and UN,t the approximation using only the first terms of the expansion. Then 



Q»[H] = E 



E 



exp(i^2u t] (fj)) 

3=1 



Error(A, /i) 
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where 



|Error(iV»| < £ E [{u tj (/,) - u N , t . 



1/2 



U = l 



The function E exp(z X)j=i u iV,i, (/?)) 
the claim follows if we can show sup n 
Lemma we see that 



is continuous in \i and Error(iV, \i) — > as N — > oo. So 
Error (A^, fi n )\ — > as N — > oo. Using the isometry as in 



E 



(u t (f)-u N , t (f)f 



H N {x,y)n{dx)n{dy) 



where 

HN{x,y) = 
is bounded by 



y - 



k=N+l 



K 2 ds 



\X} -AT 2 ! 2 



dx'dy' 



H N (x,y) < C{a,t)e- a \ x \- a \ y \{l + 



x-y 



Note that Hjq{x,y) is monotone decreasing but not continuous. The assumptions of the claim 
allow, by an approximation argument, to ignore the singularity in the function H]y(x,y) and 
replace it by a monotone decreasing continuous function Hn(x, y) of compact support. But then 
the vague convergence /i„ — > fi implies that sup^ / / H^(x,y)jj, n (dx)iJ, n (dy) j as N — > oo (for 
example by the argument of Dini's lemma). This completes the proof of the claim. I 



6.2 Density of Support 

In this subsection we give the proof of Theorem |3| ii). We start with an outline of the method. 
Assume that uo(B r (a)) > and fix T > 0. We wish to show that with probability one 
ux{B r {b)) > 0. We consider various tubes in [0, T] x R rf which connect {0} x B r (a) with 
{T} x B r (b). (By a tube we mean that for any time t the cross section of the tube with the slice 
{t} x R d is a ball of radius r.) We consider a subsolution to the equation which has Dirichlet 
boundary conditions on the edge of the tube. We will show that the probability that the subso- 
lution is non-zero at time T is a constant not depending on the tube. It is possible to construct 
an infinite family of such tubes such that each pair has very little overlap. Then a zero-one 
law will guarantee that, with probability one, at least one of the subsolutions will be non-zero. 
Applying this for a countable family of open balls we shall obtain the density of the support. 
Note this implies that the solution never dies out completely. Note also that for the equation 
(1.1) posed on a finite region the above argument fails, as there is not enough room to fit an 



infinite family of nearly disjoint tubes. 

Let us give a rigorous definition of the tubes described above. For a piecewise smooth 
function g : [0, T] — > ~R d the tube centered on g is defined as 



{(t,x) e [0,T] xH d : x e B r (g(t))} 
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If T is such a tube, let <9T be the boundary of T, minus the part of the boundary at t = and 
t = T. Now we aim find a solution (u^(dx) : < t < T) to the equation (1.1) but restricted to 
the tube T and with Dirichlet boundary conditions, that is 



— - AuJ + KujF(t, x) for (x, t) G T, 



dul 
at 

Uq(cIx) = v(dx), where supp(^) C B r (g(0)), (6.6) 

u?(dx) = for (x, t) G 8T. 

As in Section ||, a chaos expansion with respect to the noise F yields solutions to fl6.6|) . We 
do not give the proof. The only changes needed are that the stochastic integrals are restricted 
to the tube and the Green's function Gt(x — y) must be replaced by the Green's function for 
the tube G^~(x — y), that is the fundamental solution for the heat equation in the tube with 
Dirichlet boundary conditions. As in Section ^ the convergence of the series is guaranteed by the 
finiteness of an exponential Brownian bridge moments; however the moments that are needed 
are of the form 



t?i,x ,y 

0,x,y 



.2 



''" ds) /(sup \X\ - g(s)\ V \X 2 S - g(s)\ < r) 



and so are less than the corresponding moments needed to ensure the solution on the whole 
space converges. 

Fix the noise F, on its filtered probability space, and construct via chaos expansions 
{u t (dx)} the solution to ( |1 . 1[ ) started at fi G T~i a + and {u[(dx)} the solution to ( |6.6| ) started at 
v = A*lB r (g(o))- We may also construct approximating solutions u[ T ' £ \x)dx to uf(dx), by using 
the smoother noise F £ and the initial condition = Gjv, exactly as we approximated ut(dx) 
by uf\x)dx. A fairly standard comparison argument shows that u^~ ,£ \x) < uf\x). Passing 
to the limit as e — > we find that, with probability one, 

uj(dx) < u t {dx) for all < t < T. (6.7) 



We now start the proof of Theorem ^ ii). As described above it is enough to assume that 
uo(B r (a)) > 0, for some a G R d and r > 0, and to show, for fixed b G R d , that UT{B r (b)) > 
with probability one. For notational ease we shall take a = b = and r = 1 since the proof 
needs only small changes for other values of r,a,b. Let e\ be the unit vector (1,0, ... ,0). We 
consider a sequence of piecewise linear functions g n {t) given, for n G Z, by 



9n(t) 



2nte x if < t < T/2 

2n(T - t)ex if T/2 <t<T 



We write T n for the tube centered on g n . The Feynman-Kac representation (p.. 18 ), adapted for 

fT e) 

the Dirichlet boundary conditions, gives the following representation for the solution Uj, n ' (f ) 
for a test function / > supported in B%(0): 



u 



(T„,e) 
T 



(/) 



r ^°) T dx ^ e \dy)G T {x-y)f{x) 



■E, 



exp L£ F^(ds, X 8 ) J I((s,X s ) G T n , Vs < T) 
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where all the integrals are over the ball -Bi(O). By conditioning on the position of the Brownian 
bridge at time T/2 we find 



,(T„,e) 



(/) = e ~re(0)T J dx J v {e) {dy) J dzf{x) 

■G T/2 (x - nTei - z)G T/2 (nTe 1 + z - y)f(x)E 1 (y,z)E 2 (z,x) 



(6.8) 



where 



I \ rpT/2,z+nT 

Ei{y,z) = E Q ' y 



exp Uj^ /2 F( £ \ds,X s ) \ I{(s,X s ) e T n , Vs < T/2) 



and 



E 2 (z,x) =Eq£* t 



exp ( K £^F®(ds,X aHT/2) ) J I((s,X s _ {T/2) ) e T n , Vs < T) 



All the randomness in the representation ( p.8|) is contained in the Brownian bridges expectations 
E±(x,z) and E 2 (z,y). By adding a suitable linear drift to the Brownian bridge we may rewrite 



E 1 (y,z)=El / y 2 > z 



exp 



where F( n - £ )(x,t) = F( £ )(a; + nt, t) is a new noise which has the same covariance structure as 
F( £ \ This shows that the laws of E\(y,z) is independent of n, and a similar argument applies 
to E 2 (z,x) which is also independent of E\(y,z). Also for x,z £ -Bi(O) 



G T / 2 (x - nTei 



exp(— n 2 T — 2ne\ ■ (x — z)) > exp(— n 2 T — 4|n|). 



G T/2 (x - 2) 

A similar lower bound holds for Gx/ 2 {nTe\ + z — y). Using these in ( ]6.8| ) we see that the variable 



u^~ n ' £ \f) stochastically dominates the variable C(n,T)u^~°' £ \f), where C(n,T) is a strictly 
positive constant independent of e. Letting e [ we obtain the same stochastic dominance for 
the solutions driven by the singular noise F: 



,( T «), 



»/ '(f)>C(n,T)u?°\f), 

where the inequality stands for stochastic domination. Let A n be the event {u^ n (i?i(0)) > 0}. 
Then P(A n ) > P(Aq) for all n by this stochastic domination. Also P(Aq) > 0, as can be seen 
from the fact that the first moment of Uj? (Si(0)) is given, in a similar way as for the first 
moments in Ql.6|), by the heat flow in the tube and is hence non-zero. 



Finally we apply a zero-one law to conclude the result. Consider the sequence of noises 
defined by 

F n = (P(t, x + g n {t)) : < t < T, \x\ < l) for n=0,l,. . . 

Since the correlation structure of F is unchanged by piecewise linear shifts the noises {F^} are 
identically distributed and form a stationary sequence. We claim this sequence is also strong 
mixing. For this it is enough to show, for all k and bounded measurable G, H, that as n £ Z 



E [G(F_ k , . . . , F k )H(F n _ k , . . . , F n+k )\ — > E [G(F_ k , . . . , F k )] E [H(F_ k , . . . , F k )] . 



(6.9) 
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Suppose that 4>ij(x,t) are test functions supported in (0, T) x f?i(0). Suppose G and H are 
bounded continuous functions of the vector 

/ / <t>i,jdFi : -k < i < k, j = 1, . . . , k . 
JO JSi(O) / 

Each integral J / (pijdFi is a Gaussian variable. Also, the covariance between / / cpijdFi and 
/ / (fri+nj'dFi+n converges to zero as n — > oo. This implies that the mixing relation ( |6.S| ) holds 
for G, H of this special. A monotone class argument then proves the mixing relation for general 
G and H. 

Define S n to be the a- field generated by the noises (F n , F n +i, F n +%, . . .). The strong mixing 
of the sequence implies that the sigma field S = n^ =1 S n is trivial in that P(S) = or 1 for all 
S G S. The construction of the solutions by a Wiener chaos expansion shows that the solution 
U ( T «) i s measurable with respect to the sigma field generated by the noise F(t, x) for (t, x) G T n . 
Thus the event A n is S n measurable and the event {A n i.o.} is S measurable. Since P{A n ) is 
bounded below uniformly in n the event {A n i.o.} must have probability one. Finally, since 
it T (Bi(0)) > sup n n? r ' l) ( J Bi(0)) by (j6j7|) the proof is complete. 



6.3 Singularity of solutions 



In this subsection we prove the singularity assertion in Theorem |3| iii) . We first sketch a short 
argument that suggests the solutions are singular. Fix T > and x G R rf . For t G [0, T) we 
consider the process 

M t (x) = ( GT-t(x-y)ut(dy). 



It is possible to extend the martingale problem ( |1.4| ) to test functions that depend on time 
and that do not have compact support, provided that they decay faster than exponentially 
at infinity. Using the test function (t, y) — * GT~t(x — y) it follows from this extension that 
{Mt} is a nonnegative continuous local martingale for t G [0,T). The explosion principle (see 
|RW00| Corollary IV. 34.13) implies that the quadratic variation must remain bounded as t f T. 
Therefore, with probability 1, 

{M(x)) T = [ T [ [ u t (dy)u t (dz)G T ^(x - y)G T ^ t (x - z)\y - z\~ 2 < oo. (6.10) 



o 



However, a short calculation shows that if ut(y) has a continuous, strictly positive density in 
the neighborhood of (T, x) then the integral in ( |6.10 ) is infinite. 



Instead of pursuing this argument we show that the scaling relation can be used to convert 
the death of solutions at large times to the singularity of solutions at a fixed time. Applying 
the scaling Lemma with the choices a = e~ d , b = e 2 and c = e, we find that, under the 
initial condition uo(dx) = CL(dx) (where L{dx) is Lebesgue measure), that ut(B(0,e)) has the 
same distribution as e d u t / £ 2 (B(0, 1)). Also, as in the proof of Theorem |2] ii), the linearity of the 
equation and the concavity of the function z — > imply that the map fi — > Q ^[Ut(B (0, e)) 1 ^ 2 ] 
is increasing and concave in fj,. 
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Take a solution {ut(dx)} with uq of locally bounded intensity. Then, for fixed t > 0, 



E 



u t (B(0,e)) 

£ d/2 



1/2 



n 



Ql 



U t/2 (B(0,e)) 1 / 2 



?d/2 



P{ u t/2 £ dfi) (by the Markov property) 



< QE(u t/2 (dx)) 

< Qc{t)L{dx) 



1/2 



C/t/ 2 (^(0,£)) 

£ d/2 

^ /2 (B(0, e ))V2 



(by Jensen's inequality) 



rrf/2 



Qc(t)L(dx) U t / 2e 2(B(0,l)) 



1/2 



(since ito has bounded intensity) 
(by scaling) 

= (C(t)) 1/2 Q Jb(m) [C/ V2e2 (l) 1 / 2 ] (by self-duality, Proposition § 
—* as £ | (by Proposition H|). 

The same result holds true if -B(0, e) is replaced by B(x, e) for any x € H d . We may decompose 
the measure ut = uf c ^ + uf^ into its absolutely continuous and singular parts and write uf^ 
At(x)dx for a locally L 1 function At(x) > 0. Then 

i4 ac) (S(x,e)V2- 



^4y 2 (x)dx 



lim • 

e|0 



< 



< 



lim -E 



lim E 

ej.0 



£ d/2 

u^jBjx^e) 1 / 2 

£ d/2 

utiB&e) 1 / 2 



dx (Lebesgue differentiation theorem) 
dx (Fatou's lemma) 



rd/2 



dx = 0. 



Thus At = with probability one. 



In general we may decompose the initial condition no G T~C a + as a countable sum of measures 
uo = J2 n u o^ wnere each has locally bounded intensity. Use a single noise to define chaos 
expansion solutions u[ n ^ (dx) with initial conditions Uq 1 ^ . It is easy to check that J2 n u t (.dx) 
is a solutions started at uq. Then, applying the above argument to each u[ n ' yields the desired 
result in the general case. This completes the proof of Theorem iii). 



Remark: We would like to thank M. Yor for informing us of his work in [YorSC]. 
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